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ABSTRACT 


A procedure for the control system design for a small 
unmanned, untethered, underwater vehicle was specified. The 
early phases of this design procedure were performed for the 
vehicle. The linear equations of motion Hon sivall 
perturbations about an equilibrium condition were developed, 
and the values of the hydrodynamic coefficients were 
determined. After the transfer functions were developed to 
relate the vehicle motion in each of the six degrees of 
freedom to small deflections of the rudder or stern planes, 
computer programs were written to determine the poles and 


Beros of the transfer functions. The vehicle s response in 
the time domain to one degree deflections of the rudder and 
stern planes was determined and plotted. The vehicle's 


response indicated that the vehicle was stable and that the 
size of the control surfaces could be reduced to minimize 
drag unise Maintaining vehicle stability. The nonlinear 
equations representing the vehicle's motion in six degrees 
of freedom were incorporated into a computer model of the 
vehicle. The outputs of this model (¢ linear and angular 
velocities in six degrees of freedom ) are the state 
variables which will be needed during later phases of the 


control system design. modern control theory. 
Thesis supervisor: Dr. Damon E. Cummings 
Title: Visiting Lecturer 


Department of Ocean Engineering 
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Chapter I 
INTRODUCTION 
The potential uses for a small unmanned, untethered 
submersible vehicle ( SUUV) are numerous. One of the Key 


systems in the design of such a vehicle is the control 
system, which must be able to control the vehicle: 

1. during transit to and from the vehicle's operating 

location and depth; 

2. While the SUUV performs its mission on location; and 

3 during some emergencies f0r which pre-programmed 

corrective action has been specified. 

In order to perform these functions, the SUUV'sS control 
system must be able to accept pre-programmed information 
concerning course, distance and operating depths as well as 
real time inputs from various onboard sensors. The output 
from the control system is used to control the rudder and 
horizontal control surfaces and ballast control system to 
cause the SUUV to carry out its programmed assignments. 

The length of the mission which the SUUV will be able to 
perform will be dependent on the amount of energy which can 
be installed and the rate at which the available enercy is 
consumed. Therefore, the control system should be optimized 


such that the SUUV expends the smallest amount of energy 
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While it performs its mission. Designing the control system 
to minimize the overall vehicle consumption of energy 
requires that different solutions to the control system 
design problem be evaluated in terms of energy consumption. 
For example, the size of the vertical and horizontal control 
surfaces could be made large enough to provide straight line 
stability of the vehicle with no control surface deflection. 
Alternatively, the size of the control surfaces could be 
made smaller to reduce the drag and thereby reduce the 
energy required for propulsion. However, smaller control 
surfaces would require an active control system (Can energy 
drain) to provide vehicle straight line stability. 

The design of the control system should follow a logical 
sequence. Although the demarcations between different 
phases of the design process are not distinct, the following 


phases are identified for this project. 


Phase I Identity the vehicle size, shape and mission 
requirements. 
Phase II Identify a mathematical model to represent 


the motion of the vehicle in six degrees 
of treedonm. Linearize the mathematical 


model. 


Phase III Calculate the values of the hydrodynana- 


coefficients in the linearized equations. 
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Phase IV Form the transfer functions for each of the 
linearized equations. Check for vehicle 


stability by calculating the values of the 


poles and zeros for each transfer 
ENC eon. Determine the vehicle time 
response to a small control surface 
deflection. 

Phase “V Using modern control theory, model the 


vehicle in state space using the  5six 
nonlinear equations which describe its 
motion. Include in the model the effects 


of random environmental disturbances. 


Phase VI Model the outputs of the vehicle sensors 
which will be used to determine the 
GOntrol signal. Include the effects of 


random noise in the measurements. 


Phase VII Construct an observer (Kalman £1ilter) to 
estimate values of the states which cannot 
, be measured directly by the vehicle's 


SENSOEIS. 


Phase VIII Combine Phases V, VI, and. | Vil “Intoe., a 
complete model of the vehicle's motion. 
Use this model to design a linear control 
system Which will control the vehicle in 


the predicted manner. 


10 
In this thesis, Phases I through V are performed in detail, 
and a briet discussion  cH#erPhases =, VII and Vitis 
provided. 

Phases I and II are discussed in Chapter Two. Figure 1.1 
shows the principal dimensions of the SUUV. Although the 
shape and dimensions were selected arbitrarily, 1t was felt 
that these parameters would be close to those which might be 
selected for an actual vehicle design. The control system 
design procedure used in this thesis should be appropriate 
for any vehicle with missions similar to those of the SUUV. 
Figure 24 “Ls a’ -~ block diagram representation of the 
anticipated inputs and outputs for the SUUV'S control 
system. The linearized equations which are developed in 
Chapter Two are adequate to predict the vehicle's response 
to small perturbations about the assumed equilibrium 
condition of straight ahead motion with the control suriaceas 
undeflected. 

Phase III, calculation of the hydrodynamic coefficients, 
1s performed in Chapter Three. 

In Chapter Four, the linearized, small perturbation 
equations of motion for the vehicle which were developed in 


Chapter Two are solved individually to ‘find the transfer 


functions (Phase IV). These transfer functions relate 
(forward velocity), w (vertical velocity), and 6 (piten 
angle) to small deflections of the stern planes; and v 


(sideslip velocity), @ Croll angle), and w (yaw angle) to 


Jul 





FIGURE 1.1 VEHICLE SHAPE AND DIMENSIONS 
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FIGURE 1.2 CONTROL SYSTEM INPUTS ANDNOCUTFUEs 


Le 


small deflections of the rudder. The poles and zeros o¢€ 
these transfer functions were Galcwud gted. Appendix C 
contains the computer programs used to calculate the poles 
and zeros of the transfer functions. The inverse LaPlace 
transforms of the transfer functions were evaluated at a 
series of times to determine the response of the vehicle to 
a small step change in rudder or stern plane position. 
These time response curves were evaluated to determine if 
the settling time and stability of each response mode were 
Satisfactory. The computer programs used to calculate the 
time responses are listed in Appendix D. 

Having determined from the classical control theory that 
the initial estimate of the geometry of the vehicle was 
satisfactory from the aspect of stability and time response, 
in Chapter Five modern control theory was used to proceed 
with the beginning of an in-depth design of the control 
system (Phase V). Figure 1.3 shows a blocK diagram of the 
major elements of a modern control system. This thesis 
includes a development of the model used to represent the 
motion of the vehicle. This model is a set of nonlinear 
equations that describe the complete motion of the vehicle 
in six degrees of freedom —- surge, sway, heave, roll, 
Joyal 1ereai ga and yaw - and also includes terms to model the 
eftects of random environmental disturbances on the vehicle. 
A computer program is described which will evaluate the 


states (velocities and angular velocities) at a given time. 
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FIGURE 1.3 BLOCK DIAGRAM OF MODERN CONTROL SYSTEM 


We 


This program is listed in Appendix E. Follow-on work to 
this thesis should include developing computer programs to 
calculate the information discussed in phases VI and VII, 
and then using the program in Appendix E along with these 
neW programs to model the complete system ~- input, Delaticc > 
sensors, observer, disturbances, and output. This would 
allow simulation of the vehicle motion in real time and the 
vehicle response under various Conditions “could be 
evaluated. Based on these responses, the optimal feedback 
gains and observer gains can be adjusted if necessary to 
provide the required vehicle response. 

Chapter Six consists of a short summary and the 


conclusions reached by the author. 
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Chapitez in 


EQUATIONS OF HOTZON 


Zon TDENTIFTCATION OF THE VEHICLE "AND Utseeoomon 
The small unmanned, untethered, underwater vehicle for 
which this control system design procedure was developed is 


depicted aun Figure ale The principal dimensions and 


characteristics of this vehicle are: 


Length overall 11.67 feet 
Maximum diameter Zi ianches 
weight 1toz pounds 
Buoyancy 1162 pounds 


Control surfaces Rudder 
Stern Planes 
Symmetry Port-Starboara 


Deck~Keel 


The mission of the SUUV is: 

ee Transit to operating location and depth; 

Ze Operate for a predetermined period performing the 
intended mission such as monitoring the environment 
or conducting inspections of underwater equipment; 
and 


Cur Return to the predetermined location for recovery. 


a4 


Ze: DEVERORMENT “On LHe EQUATIONS OF MOTION 
The SUUV was treated as a rigid body with deck-Keel and 
port-starboard symmetry whose motion can be described using 
Six degrees of freedom - surge, sway, heave, Ba Gen.» Boi, 
and yaw. The control effectors on the SUUV include a 
rudder, horizontal control surfaces called stern planes, and 
a propeller at the stern. Mathematical equations which can 
be used to model the movement of the SUUV in all six degrees 
of freedom were then developed following the procedure used 
by Humphreys!. Appendix A contains a list of the notations 
used in these equations and throughout this thesis. These 
notations follow the standard nomenclature introduced in the 
Society of Naval Architects and Marine Engineers TMB 1-5. ¢ 


The equations which describe the forces, F, and moments, 


M, acting on the SUUV are: 


_— d tl 
eS ae (Momentum) 
VL = = (Angular momentum) (2) 


’ Humphreys,D.E.; Development of the Equations of Motion and 
Transfer Functions for Underwater Vehicles pp. 1l=15 Naval 


Coastal Systems Laboratory; Panama City, Florida, July 
197 6°. 





¢ The Society of Naval Architects and Marine Engineers, TMB 
No. i-oSreeenomenclature for Treating the Motions of a 


SUDMem@eeGesOoY Throudqh a Flusad ; Apral, 1952. 
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A set of Light hand. orthogonal coordinate axes centered at 
the center of mass, CG, of the SUUV was used to develop the 
equations of motion. Positive dizrections for the axes, 
linear velocity components, angular velocity components, 
angles, forces and moments are shown in Figure 2.1. 

The inertial forces and moments acting on the vehicle can 
be written as: 


° (3) 
2x = m(U + OW -~ RV) 


. (4) 
21 =“mtV + RUS > PH) 


° (5) 
SZ. ='wCW  tPV += CoU) 
° e (6) 
Slices (Pte = Ri + Q@R(I - I) - PEI 
x x2 = y XZ 
e (7) 
SM.-=. OT; + PRCT: “= i“) = hex + per 
y x = Xam 2 
e e (8) 
=N = RI - PI + PO(I - I ) +ERI 
Fi D6 ok yY x > BA 


The total motion of the vehiclewcanebe assumed to consist 
of two parts: Gil) an average motion representative of the 
equa librium. condrer1on,; and (€(2) a dynamic motion that 
consists of small perturbations about the average motion. 
Thus», the instantaneous velocity components at any time can 


be written: 


kg, 


SINANOW GNV SHONOA “SAILIOOTSA “SA IONV ‘SAXV JO SNOILOTNIAG AAILISOAT°?2 WNOIA 


n*y A’y 


h of ~ 
n 


@dUIIVJOY poxyy 


x 


So > = ate 








: Sree j = 
~— re 
VIUIIIJay 
TeIuoZTAoOP 


aouaIa jay 
Te uozzaoH D 


0 


x <n 
n‘y 
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Ue=" Up tau £3) 
Vo= Vo Fv Ci) 
Wo= Nowe ws C1) 
P = Po + p C42) 
2 = Qo Ff (13) 
Ric= sho az (714) 


The zero subscripts represent the average or equilibrium 
velocity components and the small letters represent the 
dynamic velocity components. By defining the equilibrium 
condition as straight ahead motion with control surfaces 
undeflected, all average velocity components except Uy are 


zero and the total velocity equations can be rewritten: 


U = Ugo + u C159) 
Vozv ( 16p) 
W=w Cie) 
P = p (18) 
Q=q (19) 


Rees taeee (20) 


ral 


Substituting these equations into the 


moment equations yields: 


Pe ont tow —rEVy) 


Sre= My +orUg + LW >. piu) 


Se NCH + pv -Guig  -au) 


2k i= pe) = or a7 (eipe Ge § 
x XZ z 
@ 
SM ee  . preCr  a-) fb) ser 
y x 2 
7 oe 
=N = ri = pr + pqaf(I -—- 
Z XZ yy. 


inertial torce and 


G2 


C22) 
4D, 
(24) 
ee) oe Cy 
y xz 
C25) 
2T ae p*I 
XZ XZ 
G26) 
a) rere (saa 
x x2 


The following assumptions will allow further simplification 


of these equations. 


es 


Assume the disturbances ‘from 

condition are small enough that 
Squares of the changes in velocity 
to the changes themselves. 

Assume the disturbance angles are 
the sines of these angles may be 


and the cosines set @€qual to one. 


the equilibrium 
the products and 


are small compared 


small enough that 


set equal to zero 


ee 


oe Products of the disturbance angles are approximately 


zero. 


The following equations result from these assumptions: 


2k = 


ait = 


2a = 


2M = 


=N = 


Additionally, using 
equations V5 20 
velocities and the 


Wei eten: 


mu 


m(v + 


m€w - 


the 


the 


rate 


(272 
(Zen) 
YUo) 
(29) 
qUo) 
. (30) 
ri 
x2 
(31) 
¥ (32) 
DE 
XZ 


small perturbation assumption and 
relationship between the angular 


of change of the angles can be 


C34) 


(35) 


25, 

In addition to the inertial forces and moments, the vehicle 
also experiences hydrodynamic forces and moments due to the 
forces exerted by the surrounding fluid. These forces and 
moments are functions of the relative velocity, acceleration 
and position as well as control surface deflections. They 
can be expressed in functional form as: 

e 8 ee e (36) 

hee 2 CU, Yi Wo vs Ws Ding ts Past G>65Vs0) 

These forces and moments can be expressed in terms o£ the 
Taylor series expansion about the equilibrium condition. 
Two assumptions allow considerable simplification of these 
expansions. 

1. Assume second order and higher terms may be neglected 
because only small perturbations from the equilibrium 
postion are allowed. 

ae Assume that since the XZ plane is a plane of 
symmetry, x ; 2 , and M are functions only of 
u,W,q,their derivatives, and @6; and 2 52 NN; and K 
are functions only of v, p>, r; their derivatives and 
p. 

Therefore, the hydrodynamic forces and moments expanded ina 


Taylor series and simplified by the above assumptions can be 


written: 
- ae spayicte + + Y RV ee Ves Sy 
Y, Y, 4s Yev fe Yv Yep + YP Y 5 3s Ss 5 (S77) 
= nr °T + +Kert+ + K,6 3 
K K + Kev + Kv + Kop + KP Kyo Ker Kr Ke (38) 


- °y +Nep + +N O+Ner+Nr+N,6 39 
Ne No + Nev + Nov Nop Nop 6? e : 5 (39) 


2h 


= + “rt a : ae 
x X Xou + Xu + Xow + xX w “+ sf + 8 + X68 + X.6 (40) 


= + ame ee nae 
Z 2, Zt 2Zut2wt Zw + a4 + Z 49 + Z,8 + 2,6 (41) 


M=MS + Meu +M ut Mew + Mw ted + Mod + M,8 + M6 (42) 


Terms similar to xX <u>%u expresses the change aneeume 
force or moment because of the disturbance velocity. Xsue 
is called a stability derivative or hydrodynamic coefficient 
and is defined as the change in the X force with respect to 


the u velocity and evaluated at the equilibrium condition. 


COX) 


u COU) 6 


The gravity and buoyancy forces and moments can be 
expanded in a similar fashion and the complete linearized, 


small perturbation equations of motion for the SUUV are: 


a Xer x -q + +X 0+ X, 6 

mus Xeut+X ut Xwt Xw + x4 x4 X, 6_°s (43) 
7 - U6) = Zu + “y 7q + + + Z. 6 | 

m(w U8) Zu Zu + 2-w + Zw ao oa o8 Z 8 8. , (ie 
rs ee a on. cae + 6 + M 6 

i Meu + Mu + Mow + Mow + ae au My — = (45) 


and 


* The notation < > indicates a subscript. 
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m(v + UoW) Taal Se ial a v6 oOR (46) 


vtY vtYeptY pty 
SN es on pP 6 


ee _ ee * ae aus in K oe K es Ker ve K r A K 6 

I ¢ TY Kev + Kv + Se a 4? 2 bp R (47) 
Sh his +NeptNptNoGOtNer+tNrtn,. 6 48 

IY 15? Nev + Nv NP af! “ye hs e: Sp R (48) 


These two groups of three equations (43-45) and (46-48) 
are functions of different variables. Hence, the motion of 
the vehicle can be examined in the vertical (XZ) om 
morreoncal CXY) plane independently of motion in the other 
plane. (The motion of the vehicle in the two planes is said 
to be uncoupled.) The first set of equations 1s usually 
called the vertical or longitudinal set of equations and the 
second set is often called the horizontal or lateral set of 


equations. 
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Chapter 22 


CALCULATION OF THE HYDRODYNAMIC COrrricie se 


S51 CONTRIBUTIONS SEROM TRESS ODy 





The values for the hydrodynamic coefficients were calculated 
by first approximating the body shape as an ellipsoid and 
calculating values for this shape, and then calculating and 
adding the contributions from the fins and the nozzle. The 


ellipsoid representing the body is depicted in Figure 3.1. 


Volume of Ellipsoid 


Surface Area of Ellipsoid 


ab 
S°= 275° 27 2s ae 
€ 


Zi 





= 5.625! 
Onc ot 


A) 
1 


eccentric Cy 


c/a 
4 


m fh oO 
T| 


(a 


Q 
i 


€ = 0.988 


PiCri oe SLLIPSOLD USED 10, REPRESENT VEHICLE 
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€ t= wer al where € = eccentricity’ 


e€ = .988 


Landweber and Johnson® determined empirical formulas to 
determine the values of many of the hydrodynamic 
coefficients. These formulas contain the terms K',K1l,and K2 
which are determined for this vehicle by using the table 
provided by Lamb®. Entering this table with the value a/b = 


5.714, the values of K', Ki, and K2 are: 


K' = 0.782 
K1 = 0.041 
Kee. =, 2052S 


Moments of Inertia 


7 eS OF for a prolate spheroid 


" Thomas,G.B.Jr.; Calculus and Analytic Geometry ; peeagoe 
Addison-Wesley Publishing Co.; Reading, Ma.; 1960. 


> Landwebe., sale and Johnson,J.L.; Prediction of Dynamic 
Stability Derivatives of an Elongated Body of Revolution ; 
NSRDC Report C—359;5 Silay, oon 


> Lamb, Sir Horace; Hydrodynamics, 6th Eda; spp. Too Dover 
Publications,New York, 1945. 
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T =m€a* + b*)7 F 
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f= 7072p r= 
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Mass 


m= Nilowy 32.2 
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m7AL 1729 LD 
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Bil 
Coefficient On Drag 
Newman’ has determined that for three dimensional bodies 
where the maximum thickness is less than one-fifth of the 
length, the — f£zactrional drag is dominant and the dzrag 
coefficient can be predicted from the flat-plate drag 
foetitacizent C<£>. Assuming a speed of five Knots, the Reynolds 


number 15s: 


Re = Ul/v- 


Re Gol 2659) ( 10209290) 7 (1.35. x) 10°) 


Re = 5.81 x 108 (for seawater at 10°C) 


From Figure 2.3 in Newman® 


C = 3.5 x 107% 
+ 


ang Since Cc = «Cc 
D Be 


Cy = 10.70035 


OX 


x = ---- 


u ou 


Ke= = 1/2 pu 


’ Newman,J.N.; Marine Hydrodynamics ; pp 20; MET Press; 
Camberdgestia.; 1977. . 
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OKT OU = == Une 


A = S§ = Surtace Area 


AR = 49.3 1* 7126.6 =s22007'- 





Z 2 = =0.2 300m! 3°] 2 - pe 


where D’ = ¢ S/l4 D =o 0 Ome 
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The circular symmetry of an ellipsoid 


causes the following coefficients to be zero: 


Sid CONTR BUPONS EROM THE FENS 








The size and shape of the fins selected was the same for the 
rudder and the stern planes. Figure 3.2 shows the 
dimensions and shape of one of these fins. The total area 
of both the port and starboard stern planes was calculated 


to be 0.4 square feet. 


The average span of the rudder and stern planes was 


calculated as tollows: 


C2578) =] 5750472 2 


WY 
i 


Sow LG 


"~ 
Tt 


The average chord was calculated to be: 


Gas OORT GY 


Aspect Ratio 


AR = S7Cc 

AR = 2.66 
Zea SOU GA. © FeO p Uo >) 
WwW £ La 
C = 27/7(01+2/7AR) 


Paso 
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FIGURE 3.2 
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nGOrmewcatcuwlating the "contributions of the fins to the 
hydrodynamic coefficients N<v>, M<w>, N<r>, N<vo; Vor>, 
wey, O.W>;, 25G>, Beal M<q>, the projected area of the nozzle 
also was included in the calculation of effective fin area 
and span. Figure 3.3 shows the effective fin used in these 


Sabewlations. 
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Area = ,.8l ft. 
Span (S) = 1.75! 
Chord (6) = 48! 





Aspect Ratio (AR) AR Svace AR =°3.65 
© 
Cr ee if 
1 + a5 


FIGURE 3.5 EFFECTIVE FIN ARBA 
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Table 3-1 lists the values of all of the hydrodynamic 


coefficients for ease of reference. 


16 


TABLE 3-1 NONDIMENSIONAL VALUES OF HYDRODYNAMIC DERIVATIVES 


Hydrodynamic 

Derivative 

BODY FINS TOTAL 

*U -0.00273 -0.90273 
oe 

u -~0.00104 -~0.00104 
*y ~0.01420 -0.025980 ~0.04018 
ue ~0.02345 -~0.000275 ~0.02373 
ue 0.00150 0.01310 0.0146 
Ye 

r 0.0 ~0.000140 -0.00014 
a ~0.01420 ~0.02598 -~0.04018 
Le -~0.02345 ~0.000275 ~0.92373 
Sa ~0.00150 Ai eies) 10 -~0.01460 
24 0.0 -0.000140 =0. 00am 
K 

=0.00302 ~0.00302 
A 0.0 ~0.0000445 -0.0000445 
Re 0 =O) C0 mol? -0.0000017 
Ne 0.01950 -~0.01311 0.006395 
M 


W 0.0 -0.000140 -0.90014 





Table 3-1 (Con't) 
Hydrodynamic 
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Chapter IV 


DEVELOPMENT OF THE TRANSFER FUNCTIONS 


Gu. 4 GENERA 
The linearised, small perturbation equations of motion were 
developed in Chapter Two. After substituting equations 33, 
34, and 35 anto the equations of motion, the longitudinal 


and lateral equations of motion can be written: 


LONGITUDINAL EQUATIONS 





mu - Xeu - Xu - Xow = Xow - ay = one = X86 = Xs °s (49) 
mw — mU 8 - Zou - Zou = Z ow = Zw - z.8 - z,° = Z.9 = 25s (50) 
- Meu - Mu - Mew - Mow - Oe - M8 - M498 = Me_°s Sa} 


LATERAL EQUATIONS’ 





av + mU - Yev = ay - use - on So yt Yo? = YY =Y, 64 (Siz) 


? R 


ll 
A 
O 


I? - Eee - Kev - RY - ae - Ss - Sky = Ket - K¥ 5 R (San 


it 
= 
O 


LY — 109 ~ Ney - Nov = Neg = Noe > N ti Nae 5, PR (54) 


The terms X¥<és>6s, Z<dsr6s, M<déds>65S, Y<dROAR, K<dR>O6R, 
and N<éR>é6R are the control forces and moments generated 


when the control system causes small deflections of the 
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stern planes or rudder. The transfer functions for the 
response of the vehicle to small deflections of the rudder 


or stern planes are developed in the following pages. 


uu .2 LONGITUDINAL TRANSFER FUNCTIONS 
Taking the LaPlace transform of the longitudinal equations 


ylelds: 


2 
=: e <= + = e = = e = = = 
{(m X-)s X Ju [- X*s - X ]w + [ zee ae X, 18 Xs °s 5) 


—_ e = = e = = e 2 = 
f Z.s Zu + [(m Z-)s Z lw + [ Zs - (Z, + mU,)s = Z.]9 = Heer (56) 


ies - eae See : . 
[ Ms M Ju + [ M-s M Jw + [(r, ays see ~ M,1¢ = M5 9s (57) 


These equations were nondimensionalized by dividing the 
force equations by 1/2 pl*uUo?, the moment equations by 172 
plt3uoZ, ena. SS. bY Uo/l. Then the nondimensionalized 


equations can be written: 





NONDIMENSIONAL LONGITUDINAL EQUATIONS OF MOTION 


2 


[ (m' = X'e)g! a Xx" ja' - [X'es' ap X' Jw! = [X'es' 
u u W W q Ss 


+ X'o'! <e = ! 
se or X,1¢ X 6 (5s) 
S 
me we! ' ' pees Lee Poets. ' 1 
[-Z 3s Zi + [(m v8 2s Ziv 


ti. 12 ' mo! ' - 7! 
5 |Z 38 + (20 +m')s' + Z,]8 = 26 °s (59) 


’ 2 ' ’ ' 
[- Mies! - Mi du’ - [M'-s' + Mi dw! + [(ty -M 2 - wee - M, 19 = M, Je (60) 
Ss 
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The transfer functions can be determined +#f060r a Wa2ven 
arnpwe by solving the transformed simultaneous equations of 
motion for the variable of interest with all other inputs 
set equal to zero. For example, the pitch attitude to stern 
plane deflection angle transfer LUMet Lone using 


determinants, can be written: 





* — wre t cyt — wt eee! ’ 
(m X ps x X 3s xX Xs 
S 
ee re ! = oa 1 _ oo? ’ 
Z “s Zh (m Z Js Zs Zs 
Ss 
- M'es' - M' = M's" ave M. 
8(s) . i y 65 
6 ,(s) (m' -_ X*-)s' = ye =X ec a x oe Xtegt? _ X's! = X!) 
u u W W q q 8 
- Z'es' - Z' (m' - Z'-)s! - 2Z' -zeg5'? - (Z' + m')s' + =z! 
U u Ww Ww q q 8 
me MiecotickMe - Mies! — mM! (I' BG o\a0- - M's! - x! 
u u W W y q q 8 
These determinants were expanded, resulting in a numerator 
polynomial ins‘ over a denominator polynomztartein s = The 


denominator polynomial for the three longitudinal transfer 
functions is a common polynomial. Setting this denominator 
equal to zero gives the characteristic equation whose roots 
are equal to the poles of the system. The damping and the 
natural frequency and the time constant of the system can be 


determined from this character2s cmc equation. The 


pul 
longitudinal transfer functions can be summarized as 


follows: 


LONGITUDINAL TRANSFER FUNCTIONS 





6 
N Z 
é ""“+Bos' + 
Wee ee OE eee e. (61) 
: avon Neon - Bs'> = cs'? + Ds' +E 
W 
Ne A ae 2 


+Bs' +Cs' +D 
Lp eS a laa ae (62) 
Ss 





‘Tong ete - Bs"? -- cst? + Ds' +E 
U 
Ni A 3"? 2 


+Bs' +Cs' +D 
se ee u u u (63) 
- Long jeu + Bs'> = cs'? + Ds' +E 





The coefficients A,B,C,D, and E and the coefficients in the 
numerators are evaluated using the equations in Appendix B. 
The values of the hydrodynamic derivatives used to evaluate 
these coefficients are as calculated in Chapter Three. iy 
order to evaluate the transfer functions, a computer progran 
Was Written which first determined the values of the 
coefficients using the equations in Appendix B, and then 
used these values to determine the roots of the numerator 
and denominator. 

Using the program VERTLIN.FORT (listed in Appendix C) the 
roots of the polynomiels in the longitudinal transfer 
functions were calculated . See Table 4Y - 1. The surge (u) 


polynomial coefficients and roots are zero because all of 


De 
OUTPUT VALUES ARE IN DIMENSIONAL FOr 
COEFFS. OF CHARACTERISTIC EQUATION 


A B C D 
O.974719E-05 0.306286E-0% 0.2062945E-04 05505 lse-c5 
ROOTS OF CHARACTERISTIC EQUATION 


oc el 
=Ui.3 041070 0.237651 
-0.361070 =0723/7651 
=0.-0 77655 0.0 


PITCH POLYNOMIAL (COEFFICIENTS 


ATHETA BTHETA CTHETA 
=0..8504S54E-05 -—0), 7547 295-05 --—0 -5s5070r—VUo 
ROOTS 
REAL IMAGINARY 
=O. 509735 0.0 
=0.,077702 0.0 


VERTICRALOVELOCETY COEFFICIENTS 


AW BW CW DW 
—-0.110338E-05 -0.463001E-05 =-0.113874E=05 —O Gy oOCe oy 


ROOTS 
REAL IMAGINARY 
=3-937068 0.0 
-0.180908 0.0 
=0.07 702, 0.0 


FORWARD SPEED POLYNOMIAL COEFFICIENTS 


AU BU CU ae 
0.0 Or 6 0.9 0.0 
ROOLS 
REAL IMAGINARY 
KKK KKK KKK KKK KKKX Ceo 
K*X ELE XE © KX MLE 0.0 
KKKKEKKKEKKKKKKKER 0.0 


TABLE 4-1 ROOTS OF VERTICAL TRANSFER FUNCTION 


E 
0.331301E-06. 
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the hydrodynamic coefficients associated with X were assumed 
to be zero except X<u> and X<u>. Figure 4.1, which is a 
flow chart for VERTLIN.FORT shows the method used _ to 
calculate the roots. Subroutine ZPOLR, a routine written by 
IMSL,Ince. and available as a library routine at C.S. Draper 
Laboratory, was used in VERTLIN.FORT to find the zeros of 
the polynomials. The pitch (@)Jand heave (w) transfer 


mUIMNCtCLORS Can bewuritten fatter cancellation): 


(S + 2.34)(S + .36 +.245)0S + .36 - 1.243) 


The vehicle's response in the vertical plane was then 
determined by finding the inverse LaPlace transform for @(s) 
and w(s) assuming aone degree deflection of the stern 
planes. Pitch angle as a function of time following a one 


degree deflection can be expressed: 


CGto m= ee O07 -— 0.0025 Cexp)(—Z. 34t) 


fo 0S (np) (30 (COs e = .065in.24t) (66) 


The vehicle's vertical velocity as a ‘function of time 
following a one degree deflection of the stern planes can be 


expressed: 


WiGt e-mecOse + .001Cexp) (C=Z2534u) 


muro OAZ(GxHD) (Se sOtICCosc. cut + .032sin.24t) (67) 
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Plots of @(t) and w(t) are shown in Figure 4.2 and 4.3 These 
plots show that the vehicle behaves ina stable, predictable 
manner in the vertical plane following a small stern plane 


wert lection. 


In the expressions for pitch angle and vertical velocity as 
functions of time, it can be seen that some terms are much 
smaller than others. This indicates that the effects of 
some poles and zeros are insignificant in comparison to the 
effects from the remaining poles and zeros. When these 
zeros and poles which have only minor effects are dropped 


trom the transfer functions, the transfer functions can be 


written: 
CS tao tl) 
Clee eS SS SS SS SS ooo (68) 
CS oo, tame 245) OS 2.36 =. 245) 
CS 18) 
WES) Stam ne een ee (69) 
CS Ome eet IOC See te Oo 2 4 ) 


The system characteristic equation can now be written: 
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Pl Vere Ole! 


9/6/1982 


PITCH ANGLE (ODEG.) 
-1.40 -1.20 -1.00 -0. 80 ~0.60 -0.40 -0.20 ~0.00 0.20 0.40 


eee 


8.00 10.00 12.00 1N.CO =—s« 16.00 = 18.00 
Tie. KSse.) 


FIGURE 4.2 PITCH MOTION FOR A ONE DEGREE STERN PLANE 
DEFLECTION 


Zev Care iar cee co Ca) 


9/6/1982 


VELOC. (FT.PER SEC) 10-3 
-1.60 -i.20  ~ -0.00 -0.40 -0.00 0.40 0.80 


Z 
-2.00 


-2.40 


“2.0% 
a oo 


3.20 


‘c.-0 2.98 ¥.09 §.c9o 12.90 18.00 16.900 18.00 20. 


8.00 10. 00 
TIME (2EE:: 


FIGURE 4.3 HEAVE MOTION FOLLOWING A ONE DEGREE STERN PLANE 
DEFLECTION 
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S2 + .76S +.187 = 0 (70) 


This equation 12S Of the -torm, 


S*¢ + 2¢wn + wn? = 0 
where ¢ = damping ratio 
wn = natural frequency 


Solving for the damping ratio, natural frequency and system 


Lime cons tanc: 


g¢ oS 0.88 
wn = Op a2 
T = Laon Secese 


4.3 DAPERAD TRANSEER FUNCTIONS 

The transfer functions for the lateral equations were 
determined and analyzed using the same procedure just 
completed for the longitudinal equations. Equations (52-54) 
were nondimensionalized in the same fashion as were the 
longitudinal equations, With the additional definition tna 


gp® = V/Uo, 


De, 


Nondimensionalized Lateral Equations of Motion 


' i t — fre as 12 fat oO t 
{(m' - Y s)s YVJB + [- Y ce ues wale 


ie Yiest* + (m' - ¥!)s"]y = ¥! 6 


aR (71) 
toa! i t ot, pe tot . yt 
{[- K oo Ki 18 + nce K ops KS Ki]¢ 
' Cyt, 7 —~ vtet - x! 
+ [(- i K Qe KS lw Ks °R DS 
fof = vy a, en \ ae (2 tot _ ay 
[- N <s Ni JB + [( ie N ae ais NA1¢ 
' t. 12 Cat - nN! 
+ [iene - N +)s NUS Jy Ns 6, C1 


R 


The lateral transfer functions were determined from these 











equations. 
Lateral Transfer Functions 
nf 3 2 
6 s'(A.s'~ + B.s'” + C.s' + D,) 
PRS B B B B 
B/S, 7 7 4 3 2 C74) 
atae s'(As'’ + Bs'~ + Cs'” + Ds' + E) 
> 
ae s'(A ar +B s' +C,) 
o/b, = Se eS a 75) 
x ‘Lat Sugeest - Bs'> + cs'* + Ds' + E) 
y 
i Ast? +B s+ s'+D 
~/6. = ee el eee (76) 
S oree Sr nee + Bs'> + or + Ds' + E) 
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Figure 4-4 is aflow diagram of the computer program 
HORLIN.FORT (listed in Appenaz =e e The values of the 
coefficients and roots for the lateral transfer functions 
were calculated with HORLIN.FORT and are listed in Table 
U—-2. The coefficients and the roots of the roll polynomial 
are equal to zero because of the vehicle's symmetry in both 
the vertical and horizontal planes. The transfer functions 
for the sway and yaw motion are written (after 
cancellation): 


(S + 4.38) 
Vis)/6xr(S) = eee eee ------------- (77) 


WiCS TAO CS). 4a i8 a cee oe ee ee (78) 
CS Fe ORS: ed) 


v(t) and w(t) were determined: 


~ 


v(t) = 0.58 - 0.62exp(>. 45 (79) 
v(t) = 1.014 —= 0.010 > 0.0 Zero (80) 


PLOeS: “Of. VCE) and w(t) are shown in Figures 4.5 and 4.6. 
These plots show that the vehicle is stable and behaves in a 
predictable manner in the horizontal plane tor ‘smaill 
deflections of the rudder. 

In order to maximize the use of the installed energy, the 
control surfaces should be redesigned to reduce their area, 


which would reduce the drag on the vehicle and thereby 
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OUTPUT YALUES ARE IN DIMENSIONAL FORM 
COEFFS. OF CHARACTERISII©G EQUATION 


A B Cc D E 
0.301377E-07 O.1TOS5S44HE-06 O0.899080E-06 0.257019E-05 0.9765227—05 
ROOTS OF CHARACTERISTIC EQUATION 


-0.264699 5.186463. 
-0.264699 -5.186463 
"4-07235c6 0.0 
-0.446244 0.0 


SIDESLEP POLYNROMDARE VeOrRFES. 


BETA BBETA CBETA DBETA 
0.279417E-07 0.137099E-06 0.818321E-06 0.329856E-05 
ROOTS 
REAL IMAGINARY 
~4 377226 0.0 
-0.264695 =-5.186463 
=O 264695 5.186463 
ROLL POLYNOMIAL COEFFS. 
APHI BPHI CPHI 
0.0 0.0 
ROOTS 
REAL IMAGINARY 
¥¥ YE YK KYM EK KE KK YX 0.0 
KKK KKK KKK KKK KX 0.0 


YAW POLYNOMIAL COEFFS. 


APSE Crom Pe owe 
0.2 20054 50n—0 = Oy 5.002 -0.738029E—-06 > —-05 579 1 SGUEgucG 


I 
E07 
Regis 
RPEAL IMAGINARY 
=O). 2.6 46.97 Seo 465 
-0.264697 =5 186465 
-0.814914 0.0 


TABLE 4-2 ROOTS OF HORIZONTAL TRANSFER FUNCTIONS 
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OWAY MOTION 


Shea)! JUS es 


5.40 6.00 


4.80 


(ETSPER 


SWAY VELOC. 
1.80 2.40 3.00 


-20 


0.60 


8.00 10.00 12.00 14.00 16.00 18.00 20.00 
ly ae Glee ae 


EOeGURE 425 SWAY MOTION FOR A ONE DEGREE 
RUDDER DEFLECTION 


YAW MOTION 


3/6/1S82 


YAW ANGLE (DEG.) 
-0.00 -7.20 -5.60 -4.00 “2.40 -0.80 0.80 2.40 4.00 5.60 


ie ooo oe 


~10,40 


ca aca Ses E-ce o 10.39 12.00 14.00 18.00 = 18.09 2c 
"sh Fay 


FIGURE 4.6 YAW MOTION FOR A ONE DEGREE RUDDER DESLECria 
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maximize mission range or time. Atter selecting a new L£1n 


size and shape, the affected calculations in Chapter III 


should be repeated oO determine new hydrodynamic 
coefficients. Dene ne wana lyctc Ot matnacmchapter should be 
repeated. This sequence should be repeated until the size 


of the control surfaces is the smallest which will provide 


the desired control characteristics. 
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Chapter V 


APPLICATION OF MODERN CONTRO REOR: 


After analyZing the vehicle uSing clasSical control Jtheor 
the next step was to apply the principals of modern control 
theory which uses state space, state variables and state 
vectors to describe the dynamic system which is to be 
controlled: Modern control theory can be used to design 
control systems which have multiple inputs and multiple 
outputs. 

The state variables are a set of variables of interest, 
X<i>, which completely describe the state of the system at 
_ any fixed time. The finite number, n> of state variables 
which can completely describe the system at any instant form 
an n component state vector: 

x =. CX 27 hoes ASOD 

The inputs to the system are represented by a control 
vector, U. 

The state equations which represent the system are then 
Written in the form: 


at 
x 
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ne Lem Ahe co sama Dp are coOetrtICilent matrices and YCt) is the 


Sutput vector . The state variables are not aunique set. 
That is, various sets of state variables can be used. Where 
possible, it 1s advantageous to choose as the state 


variables those variables which have physical significance 
and which can be measured. 

The next step in the design sequence was to model the 
nonlinear motion o£ the vehicle in six degrees of freedom. 
Appendix D contains a listing of the six nonlinear equations 
Which were proposed as the standard equations of motion for 
submarine simulation®. 

Append1x @ GoueRbne the listing o£ a computer subroutine 
which uses these equations to calculate the values of the 
vehicle's linear and angular velocities in both the body 
coordinate system and the earth coordinate system. The flow 
enart £or this subroutine is shown in Figure 5.1. 

This completes Phases I through V of the control system 
design process as discussed in Chapter I. Addataonak, work 
on this topic should begin with Phase VI and procede to the 


einal control system design. 


Gerter, M. and Hagen, G.R.; Standard Equations of Motion 
for Submarine Simulation; NSRDC Report No. 25.0% 
Washington D.C. ; June, 1967 


10 


Lee, Jang Gyu; The Charles Stark Draper Laboratory, Inc.; February, 1982 
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Chapter VI 


CONCLUSIONS 


A general procedure for the control system design of a small 
unmanned, untethered, underwater vehicle was formulated. 
The early phases of the procedure were accomplished for the 
vehicle. accomplished for a vehicle. 

The linearized equations of motion in six degrees of 
freedom for small perturbations about an equilibrium 
condition of straight ahead motion were formulated, and 
values for each of the hydrodynamic coefficients in these 
equations were calculated. The vehicle's response to a 
small control surface deflection was analyzed by forming the 
transfer functions f0r each of the equations and then 
solving for the vehicle response in the time domain. The 
initial design of the control surfaces for the selected 
vehicle size and shape was conservative. Smaller sconcro. 
surfaces should be designed and the analysis should be 
repeated until the smallest control surfaces which provide 
Satisfactory vehicle control are achieved. 

A model of the vehicle was formulate using the nonlinear 
equations of motion. Additional work in the area of this 
thesis would be to complete a nonlinear model of the vehicle 


and its control system. Modeling of the vehicle's sensors 
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and measurement noise terms and constructing an observer to 
calculate non-observable states are required to combine with 
the nonlinear vehicle model and input terms. 

AU) gts likes complete nonlinear model should then be used to 
predict the vehicle's response in real time. An actual 
control system should then be designed that has the desired 
characteristics to control the vehicle within the specitied 


tolerances. 


Symbol 


CB. 
EG 


I 


yz 


K 
Plp| 


Pg 


pal 


APPENDIX A 


NOTATION 


Dimensionless Form 


B 
pt*ue 


wo 
| 


ne 








Definition 
Buoyancy force, positive upward 


Center of buoyancy of submarine 


Center of mass of submarine 


Moment of inertia of submarine about x axis 
Moment of inertia of submarine about y axis 
Moment of inertia of submarine about z axis 
Ercduct of inertia about xy axis 
Product of inertia about yz axes 
Product of inertia about zx axes 


Hydrodynamic moment component about x 
axis (rolling moment ) 


Rolling moment when body angle (a, 8) and 
control surface angles are zero 


Coefficient used in representing Kz as a 
function of (n-1) 


First order coefficient used in gad at 
K as a function of p 


Coefficient used in representing K as a function 
of’ p 


Second order coefficient used in representing 
K as a function of p 


Coefficient used in representing K as a function 
of the product pq 


qr 


K 
Vv 


“*'p 
wr 


ér 


M 
47 


M 














2 


Coefficient used in representing Kas a 
function of the product qr 


First order coefficient used in representing 
K as a function of r : 


Coefficient used in representing K as a 
function of r 


First order coefficient used in repréesentime 
K as a function of v 


Coefficient used in representing K asa 
function of v 


t 


Second order coefficient used in representing 
K as a function of v 


Coefficient used in representing K as a function 
of the product vq 


Coefficient used in representing K as a function 
of the product vw 


Coefficient used in representing K as a function 


“of the product wp 


Coefficient used in representing K as a function 
of the product wr 


First order coefficient used in representing 
K as a function of 6, 


Overall length of submarine 


Mass of submarine, including water in free- 
flooding spaces 


Hydrodynamic moment component about y axis 
(pitching moment) : 


Pitching moment when body angles (qa, B) and 
control surface angles are zcro 


Second order cocfficient used in representing 
M as a function of p. First order coefficient ts 
zero. 


First order coefficient used in representing 
M as a function of q 


First order coefficient used in representing ~ 
sil as a function of (n-1]) 


Coefficient used in representing M as a 
function of q 


M 
ala] 


Mi alos 


< 
" 
z 














HE 


Second order coefficient used in representing 
M as a function of q | 


Coefficient used in representing M,_ asa 
function q a 


Coefficient used in representing Masa 
function of the product rp 


Second order coefficient used in representing 
M as a function of r. First order cocfficient 
is zero 


Coefficient used in representing Masa 
function of the product vp 


Coefficient used in representing Masa 
function of the product vr 


Second order coefficient used in representing 
M as a function of v 


First order coefficient used in representing 
M as a function of w 


First order coefficient used in répresenting 


M., as a function of (n-1) 


Coefficient used in representing M as a function 
of w 


First order coefficient used in representing M 
as a function of w; equal to zero for symmetrical 
function 


Coefficient used in representing Mg as a function 
of w 


Second order coefficient used in representing 
M as a function of w 


> 


First order coefficient used in representing 


w| w as a function of (7n-1}) 


Second order coefficient used in representing 
Mas a function of w; equal to zero for sym- 
metrical function 


First order coefficient used in representing 
M as a function of by 


First order coefficient used in representing 
M as a function of §, 


First order coefficient used in representing 
Msg. as a function of (7-1) 


74, 








N N' = Bot uz Hydrodynamic moment component about z 
Be axis (yawing moment) 
Nee 
Ny Ny = ZU Yawing moment when body angles (q, 8) and 
p control surface angles are zero 
N N! mer F ff 3 N 
= irst order coefficient used in representing 
P P Pe as a function of p 
N. N.' = ~t Coefficient used in representing N as a function 
Pp Pp pt of p 
‘pq 
N No = ra Coefficient used in representing N as a function 
Pd Pa tpt of the product pq 
N Reese Coefficient used i ting N functi 
qr qr = ple oellicient used in representing as a tunction 
of the product qr 
N 
N N_'= — First order coefficient used in representing N 
r r $pt4U 
as a function of r 
Ny au : 
N N ‘= tot First order coefficient used in representing 
a iu peru Ny as a function of (n-1) 
N: net : 
N. N.' = —t— Coefficient used in representing Nas a function 
D9 r s £& ° 
p of r 
Nr{r| ) 
N N ie Second order coefficient used in representin 
yy “ae Pp & 
=F rir] ~ gpt N as a function of r 
N 
Niel br Niet a, = 2 BE Coefficient used in representing N;.. asa 
$pL4U function of r 
Ny . e¢ . ° 
Ny Se 3 First order coefficient used in representing N 
ep" U as a function of. v 
N 
N N '= -—L First order coefficient used in representing N, 
ac pa pt” U as a function of (n-1) 
Ny _ , . 
N. N.'=—— Coefficient used in representing Nasa 
Vv Vv &oL* . e 
function of v 
Ny : ; 
Nn Nua ——t Coefficient used in representing Nas a function 
a q spt‘ of the product vq 
N 
Nii a Male = as Coefficient used in representing N; asa 
tot function of v 
N 
ae ay al = slut Second order coefficient used in representing 
Spt N as a function of v 
INern et 
Nvlvin Ne) vin = ee First order coefficient used in representing 


el as a function of (7-1) 


érn 


tte 


Qe 


be 


ce 


de 
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eile 


ea 
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Gale ce 
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Coefficient used in representing Nas a function 
of the product vw : 


Coefficient used in representing N as a function 
of the product wp 


Coefficient used in representing N as a function 
of the product wr 


First order coefficient.used in representing N 
as a function of $y 


First order coefficient used in representing 
Ng, as a function ef (n-1)- 


Angular velocity component about y axis 
relative to fluid (roll) 


Anpular acceleration component about x axis 
relative to fluid 


Angular velocity component about y axis relative 
to fluid (pitch) 


Angular acceleration component about y axis 
relative to fluid 


Angular velocity component about 2 axis 
relative to fluid (yaw) 


Angular acceleration component about z axis 
relative to fluid 


Linear velocity of origin of body axes relative 
to fluid 


Component of U in direction of the x axis 


Time rate of change of u in direction of the 
x axis 


Command speed: steady value of ahead speed 


. componcnt u for a given propeller rpm when 


body angles (a, 8) and control surface angles 
are zero. Sign changes with propeller reversal 


Component of U in direction of the y axis 


Time rate of change of v in direction of the 
y axis 


*B 


=G) 


Xx 
aq 


rp 


rr 
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Component of U indirection of the z axis 


Time rate of change of w in direction of the 


Zz axis 


Weight, including water in free flooding spaces 


Longitudinal body axis; also the coordinate ofa 
point relative to the origin of body axes 


The x coordinate of CB 


The x coordinate of CG 


A coordinate of the displacement of CG relative 
to the origin of a set of fixed axes 


Hydrodynamic force component along x axis 
(longitudinal, or axial, force) 


Second order coefficient used in representing 


X as a function of gq. 


is zero 


First order coefficient 


Coefficient used in representing X as a function 


of the product rp 


Second order coefficient used in representing 


X as a function of r. 


zero 


First order coefficient is 


Coefficient used in representing X as a function 


of u 


Second order coefficient used in representing 
X as a function of uin the non-propelled case. 
First order coefficient is zero 


Coefficient used in representing X as a function 


of the product vr 


Second order cocffient used in representing X 


as a function of v, First order coefficient is zero 


First order coefficient used in representingaaa 


as a function of (n-1) 


Coefficient used in representing X as a function 


of the product wq 


Rebs 
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ee 
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t Y 
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eee = 2 
gpt°U 
ty 
Y's 
P €pl°U 
eee 
es 
x '= aD, 
p|p! spt 


a 


Second order coefficient used in representing 
X asa function of w. First order coefficient is 
zero 


Pirstworder coeiiicientarsed in répresenting X,,., 
as a function of (7-1) 


Second order coefficient used in representing X 
as a function of 6). First order coefficient 
1s Zero 


Second order coefficient used in representing 


Mas a tunection of oe First order coefficient is 
zero 


First order coefficient used in representing 


X br asa function of (7-1) 


Second order coefficient used in representing X 
as a function of 6,. First order coefficient is 
zero 


Poirscmonrder scoelticient Gsed In representine 


Xe sb as a function of (7-1) 


Lateral body axis; also the coordinate of a 
point relative to the origin of body axes 


The y coordinate of CB 


The y coordinate of CG 


A coordinate of the displacement of CG relative 
to the origin of a set of fixed axes 


Hydrodynamic force component along y axis 
(lateral force) 


Lateral force when body angles (a6) and control 
surface angles are zero 


First order coefficient used in representing 
Y as a function of p 


Coefficient used in representing Y as a function 
of p 


Second order coefficient used in representing 
Y asea function of p 
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Coefficient used in representing Y as a function 
of the product pq | 


Coefficient used in representing Y as a function 
of the product qr 


First order coefficient used in representingss 
as a function of r 


First order coeificivent used in "represen 
Y, 2s a function of (74-1) 


Coefficient used in representing Y as a function 
of r 


Coefficient used in representing Y,_ asa 


function of r 


ér 


First order coefficient used in representing 
Y as a function of v 


First order coefiicient used in represemiine 
Y, as a function of (7-1) 


Coefficient used in representing Yasa 
function of ¥ 


Coefficient used in representing Y as a function 
of the product vq 


Coefficient used in representing Y, as a function 
OL x 


Second order coefficient used in representing 
Yas a function of v 


First order coefficient used in representing 
Yviv] asa function of (7-1) 


Coefficient used in repre enting Yeacea 
function of the product w 


Coefficient used in representing Yasa 
function of the product wp 


Coefficient uscd tn représcnting Yas a 
function of the product wr 


First order coetficient used in represenem. 
Y as a function of ér 


First order cocfficicnt us@d in representing 


os 
Ys, as a function of (n- 1) 
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Normial body axis; also the coordinate of a 
point relative to the origin of body axes 


The z coordinate of CB 


The z coordinate of CG 


A coordinate of the displacement of CG 
relative to the origin of a set of fixed axes 


Hydrodynamic force component along z 
axis (normal force) 


Normal force when body angles (a, B) and 
control surface angles are zero 


Second order coefficient used in representing 
ZeaGea 1UnGiiOn+ 1p. First order cocificient 
is zcro 


Paxrst order coefficient used in representing 
Zasa function of q 


bars: order cociiicient Used in representing 
Z, as a function of (7-1) 


a 


Coefficient used in representing Z asa 
function of q 


Coefficient used in representing Be 2s a 
humelIOn Os q 


Coefficient used in representing Z asa 
function of the product rp 


Second order coefficient used in representing 


Zasa function of r. First order coefficient 
is zero - 


First order cocfficient used in representing 
Z as a function of w 


First order cocfficient used in representing 


Zu, as a function of (7-1) 


Cocfficient used in representing Z asa 
function of w 


First order coefficient used in representing 
4 as a function of w; equal to zero for sym- 
metrical function 


Coefficient used in representing Z asa 
function of q bs 


z ? 
2 7, 2 w [| w | 
wiw] wi w| pL? 
7 
_ owl 
Zuatwin Zwiwin > Eps? 
Zz 
WW 
Zw Zw * Epes 
1b, 
rip (ogee Sie 
Z Z ‘ts a abs 
§s 6s $pt*u? 
Zz Zz. +2 78s 
bs 6s) pt ?U? 
a 
B 
*, 
Ses 
$s 
7 
8 
? 
= no ee = 
la ne 
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Second order coefficient used in representing 
2 asa function of w 


First order coefficient used in representing 


bai cal as a function of (7-1) 


Second order cocfficient used in representing 
Zasa function of w; equal to zero for sym- 
metrical function 


First order coefficient used in representing Z 
as a function of 6, 


First order coefficient used in representing 
Z as a function of 6, 


First order coefficient used in represencme 
Z,_ as a function of (7-1) 


§s 
Angle of attack 
Angle of drift 
Deflection ef bowplane or sailplane 
Deflection of rudder 
Deflection of sternplane 
a. =e 
The ratio —— 
U 
Angle of pitch 
Angle of yaw 


Angle of roll 


Sets of constants used in the representation of 
propeller thrust in the axtal equation 
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APPENDIX B 
EXPRESSIONS FOR THE LONGITUDINAL TRANSFER FUNCTION COEFFICIENTS 


The longitudinal characteristic equation is 


= ete + Bs'> + cst? + Ds' +E 


sone 
where 
A = (n' - x") (a! - at Gy = fe - X's aes eee = Zt Ms oe 
- (n' - A Ne - ee BO - ee - M's « (m' - x24 : 
B= - (n' - X's) Cm! - a - Z 7 (m! = mY Gy -~ we) 


- X'(m' - Z'+)(I' - M'-) = X'oM'e a @q +m') - Mi X'e mas - M'eX'Z'- 
u wy q Ww q u 


w eq 
— 7remMtey! — 7IemMiot. — 7Iinteyt. — ' _ ' ' ' ec 
Z ead Z 3M 2 q ZM ok 4 (m "2M fog. Mi (am Z =) q 
+M'eZ'X's + Z'eX' om! - ZX (I' - iy - Z'*X'(I' - M's) 
uw q u owq u wy uw sy q 


Praag? ee, a ' ‘) ~ mt ' cyt, foe mrex?Zte . 
M ¢ (m X a (24 + m') Mi (m X 2 q M oxy 4 


a ' _ yt, t 9, ' t/t 2 yt, ' tft — 7t, t 
C (m X “(a Z Mg = Z/(m X fe =< Xi (a z ys 


' oOo mt.) o yt emte7! < Miyte(7! ') ~ wiexyt(7? ' 
ot 2A, M a X aM 526 M4 oe 710) M bw \2q + m') 


— Z'oM'+X! — Z'eM'X!' = Z'M'eX! = Z'M'X's = (m' = Z' «)M' +X 
uw 6 uw q weg uw q u 6 


= t eee ? ? . 4 a ee ame (oe “Xx'm' 
M (a Ze a eel eo aa Le m8 We q De ee wig Zz. we w'q a weg 


= fee rhe Ns agli fet 2 wee) iat ay 8, Jee 
fe ae M ay M « (m X 26 M) (m X i? m’) 


jee ex’ (Z' +m’) + X'M'Z's = M'X'Z*e . 
wu’ q uw q uw q 


D= 
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Mi (m! — xte) 4 ¥' fm’ = Z'-)M! = ZINIM! | mist ee eo 
2M, (m x a Xi (a Z Mg aoe Mix 22 M oX2 





@ 8 
—~ Miyt(¢7? ') © 7temiy!? © gimtoyt — 7iwty! . wi¢m!t — 7t.\y! 
sre + m') z oM%e ZoM oe A Pit M (m Z. 3X 
Le i] ' + t L t+ L] tome row t ? i] Vey es ¢ ore eS ' 
+M ewe Aes ZX Me + Z ate + ee M(m X “26 
+ M'ex'zZ' + M'X'(Zz' + m') 3 
wou @ wu q 
Erm - Z2'X'MS. = M'X'Z! = Z'M’X! + Z'M'X! 4+ Ze 
wu @ uw 6 uw 6 wu 6 wu @ uw @ 
The pitch response transfer function is 
Q 
5 A ue + BS ae 
Ss 8 8 8 
1/5. = = 
Long Long 
where 
= ' ms ie ve tig ’ tewte a ye t'eomt. ' UP, oh 
Ay oes X “cm Z ) + ase aM = *5 7 M ‘gl none 37 . 


' t _ at, t. ' fe ' yt, 
+ X' (m Z eM a +Z'M ¢ (a Xx a 


Co 


6 


6 
e e 


_ mi ' _ ytey7t o wt ytf(mt _ 7. ' yt.wt t yin. 
oe X ae ea Z ne Bs on oN + “6 . 


+X' Z'eM' + x! Z'M's = M! x'eZ' = M! X'Z'> + X! (m' - Z'?)M' 
u W u Ww wou w ou oa wou 


6 6 6 6 
e e e e 


' 7twt, - ot wtyyt MY fer oe Pee ae 
3 ow sone ony + 75 ws X 7 ; 


=M' x'Z' + Z' X'M' + X! Z'M’ = M! X°2" = X77 ee 
u Ww wou uw wou wou a. wou 


6 6 6 6 6 
e e e e e 


The vertical velocity transfer function is 
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Ww 
; ae A ale + B aie +Cs' +D 
Wwe S - Ww Ww Ww 
95 ‘Long Bione 
where 
= ' ree t. oe: 1. ’ te.7t, t PS ae ’ ert ae 
A, Zs (m X wy M oe 3 X« M 52 , as M. Z SX q Zs M ik q 
e e e e 
+ ] Ue ef) = Uo t i. ur Uc : 
a ey M aU a5 eg X 2 : 
=> - ’ = t, ee ' ' t= t. ' tot, ' i. joe ' 
a Zs (m X ee Zs 8 Oa. M =) — Xs MZ e ot X M ee nm’ ) 
e e e e 
=e Ms Zh +X! + Ms Zikl: - Zs Mi ex! = Zs Mix’: - X. Z' MT 
5 q % q a q a gq 5 gq 
? ’ eens i. ' ' o t, ? ' Zz ' tot, 
+ X; Ae ai M 2 + Ms (m X we ee + m') Ms X12 a 
e e e 
ran ' ' t. ' ' tut ’ tot ' ' ' ' ' teyt 
C., Zs (m X Mg + Z 5 ee + Xs M ° 29 + Xs ye +m') + Ms Z ee 
e e e e e 
' tot ' tevt — ' twtr ’ tem? To ' tut ' ee ue t 
ae M. fae Zs M 58 Zs ina X« Z ig X bea + M. (m X 26 
e e e e e e 
_ ' ' ’ ' 
Me ee +m'). 
e 
Ny 2h tut ' i ' ty! - ot i a a 'wt om! tot 
te Zs X Mg os X M 29 + Ms Z xe Zs M io X. Z Mo Ms X 26 
e e e e e e 


The forward speed transfer function is 


: ae A 5? + B ate +cCs' +D 
u u u u 


A 
S ovens Long 
where 


= i] t= 5 oe i. ' f.5t. ? tout < ' ? 
A, = Xs (am BN Ne Mane ine cogs got is 


S e e e 


' t. f = w = ' '.7', 
+ Aye Ste M 2 aoe ae G 


SARS DOr 
we q 


or 


| t ti oe oe t Ll. Cy ' ' t. ' ' ' 
X, (m Z wg Xs = M alas + Ma x ae +m') + M. 
c 2 € s 
+ Yb M'ex'! + Y died M'x'te + M' (m' a Lok = M! Vf GET mae Vie X' em! 
6, wdq Sow 4q 6. wq ced 6, w4q 
+Z' X'(1' - M's) - Xi M's (zZ' +m") - Xi Mizts . 
os y 20 oe a q ) 6 w q 
a is ' one te ’ ’ EVA) ’ ! 7! ' ’ ' ? 
= x, (m Z Mo + Xs ZG +e Ms X “26 oa M. eS +m') + Z 
e e e e 
’ tyt ' a ae i om ' ' a ’ aes ’ Trt 
+ Zs HG + Me (1 Z Xp Me Z aa Zs Xx wg Zs aa 
e e e e < 
_ x? '.7!' _ yt ' '.2. np! 
oe a Ki Ma m') 
= X' Z'M' +M! X'Z' + Z° MX? —- Ml ZX Sea ee 
w 6 Ww w w 0 a w w 6 


6 
e 


6 8 6 8 6 
e = e 


8 6 
e 


» lee 
wW 


6 
e 


q 


q M'*X 
WwW 


8 
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EXPRESSIONS FOR THE LATERAL TRANSFER FUNCTION COEFFICIENTS 


The lateral characteristic equation is 


A = BOO + Bs'> 1 


+ 
Lat es 


+ Ds' + E) 
where 


= f 2 x7 te. fa aie 1 eee he < ee, weer it, ie, 
A (m .) (1, N sHCE K y. aoa a I. K 2) b 


+ Y'eK'e(- 7 ' = N'e) = Neve (Cl = K'e) ee K'-(L' = N'e)Y'e 
Ey KZ P Ven x Pp Tame, rp 


Pe Ty = WT, ~ FD 
Be = Ont = Wel TL SAMUS GD) — ee eR — |:o) 


- Y'(I' - N'-)(E' - K'+) + N'e(- T' - K's)y! - NheKty'- 
WZ Coe. x P Vv XZ r’ p vr p 


“+ 


N'(- i = Ky) ys Pe Y'eK' en! + Y'-K'(- Th = N's) 
Vv oz rp rovp raw XZ >) 


as ve ' ieoa eee t = i; wre v% q uae eo ? : am eae ; 
(m ap epee ) pee enraaaN + (m ae = K » 


tvt, feet toy ete er ay ts t tontver. o t _ wt. tvwt, 
NY er K 2 K o(1, N as saa oN LY 5 (I N KY 2 


!' _ vit.yfo- Tt _ xt. ' ! wt, ‘7. Tt cw, 
+ (m yé oi Ife K eke + (m Y )K I N 


ie aa ¢ = UE = ' - '. 
e i ie : aK tz y = 


The 
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t i te 2 eee a i a fe ' ' 2 Ve ’ 
Nm i DES (m x ely N Kg + Yi (r N os 

vi f. wte eg ae t . yt. 1 ~ w_tervive Oy A on et) 5 ee ' 
NAY CL, K - + N aK Lae K Y 4 N vroe + NYC _ K Bhat 


N'K'Y'= = YSN) = Y"eK'NY 47 Gn! =) NE 
vr Dp rv? rv p r "ve 


(m' - Y')K'(- I'. - N's) + No*Y'eK! = N's Gi eK ene 
st, XZ P VE 9 Vv ip Vga 


‘r,t _ vi '" ~ vtey ._yvte/fT? — wte\vi tentivir _ ' ~ nte\ytyt 
No (m YW), K y. K (1, N oY ork Pes (I. N 7 


tn tv. 4 oe SPO YN GS OP ee ee Nar ' ~ vtleyvin’? 
KEY A (m B 4 va) ies K DN, (m Be ae 


Y'(- tT - K'-)N' - K'Y'(- L* - N's) 
Vv XZ r Pp a XZ Pp 


N'(n' eed y's)! of. Yoci- a N'*)K! ae Neve i — NieKkty' 
6 vi voz ro? bev op vr 9 


ie TU Ses yy ULoONely 1 oe. y tent ' _ vi\yvtent 
shal Ie K oY, eee E ve Ms (m YK we 


fo Vvtyetnt = Nle(m! ~ Y'yr' + n'y'exw® -N' fp | YR) ee ene 
(m SE N <(m YK, NYY aS N (a eS K yy 6 





 -— nt tw! tnivi — ' _~ wte tn! 2 vif rt o ye ' 
(I) N WKY, + SEL (m v6 pkg bes tx2 K NG 
K'y'N’ 
rvp. 
—~ Nivirt —whirv'tyv't foy'yr'n' = N' tm! = VY")! + KIN’! Se 
See BE poet. + (m aS N(m eS eae KY UN 
sideslip transfer function is 
Nj "(A a'> + Bo) cee pena) 
ie tee B B B B 
SR alae Tate 


ii 
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where 
= ’ [' = ve eee '. = ' a ' oe i. gee ’ Pee (pee ' Ey i 
A, GA eS oe. K +) ws I’. K yy a coe Y aC I N'e) 
+ ys Ve = i: ’ a te pee ’ a ’ = . = ’ = '. 
Ne AGUs K +) + K; (1) N a : Y3 ( i K aK sl N a 
Ir Ir Tt 
a ’ = fs ‘ ’ ' es %. =m ’ _ ’ se Fy ’ ’ tw, 
Be Y; (I! N WK, Ys Ne K y Ne ( eeu K ie + Ni KIY 
rt rT aS 1 i 
se ' - + ’ = ’ i ’ = ', = ’ A an 
Ks | av Ky (m' reo i N *) NS Y KY 
rt 22 
= (m’ = Y' Ng mee - Kt. SY + Ks (1% - Nfe Yo - K'N ah. 
6 Pp 
= Ir 
’ a ’ = i ? ’ = ’ = '. 
+ ws ( ne K iN) + ¥, K aK oe N 2 
Tr Ir 
a ’ eas en ’ ’ ’ a ’ a 2 ’ ’ 
Ce Y; gt, N Dy + YS N! kp oe a K wYy ~- Ne , oe 


ook. Y°sN? = Ks (m' - Y')N' - NE 
r Pp 


‘oy a ' 
5 Ny Y'eK' + (m' ye DNs K 
‘ oF 


= r , P 


: ' 
5 K! ro . 


T 


a ee eae Sa ' ' ' ‘(2 Tt 2 te ' 
- Bea : N Se - N! oss -P a Ys ., I. K DN, 24 


= ’ i] ® =) ys *\nr7? = trem! — te wt 
Da Soa mo + N« K! rg Ke aa YONG + (n' » Ns 4 NUK 4 


- o KIN, 
r 9 


The roll transfer function is 





) 
N De 
bp s'(A,s'" + B,s' + C,) 
> . pe et Ou 
SR he AS oe: 


where 
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« ' ee t, 17 Hos - ' ey ’ = Veo ’ vow th 
i Mo vs ss ee os : oy we : a ae 3s ey Bae s ap r 
r r & 
= ’ ee, © 2 a ’ he a a = ' te '. = ' = oe 
Ke Ni ey's + YY, KC GAN 2) ee Y Gat, 7 ee 
| Gi r ps 
pane ’ ee toes ' ' oes ts ' 'svt ’ Vi fo ' = ie 
ay Ks (m Ms =) Ne K. Y(t, N Ey + Ys N a Ys No le K 2) 
: 1g r r 
= ' V4 pl heas ’ ’ as a ’ ie 4 = ' as ' Vy (ee ' teornrt 
ce + (m Y ou Ns Ave = 7 ‘oo + (m ie) “svt 2 ne yh 


+ 


' ' —_ ', ' se te ’ ’ of oe 0) ere 
Y Kl, N ) +0N 20m a6 SKS +i Yat ie K =) 


6 6 6 
r r r 
- xv' yviwt 'wiyvt . wnt vwtf(nt — vi 'wt¢at . vty . vt yt 
a Ks YN, 5 Ys NRE No Ky (mn ue) + Ks N(m te) Ys KN, 
r r r 1g rE 
~ wt vir 
N« YK. . 
r 


The yaw transfer function is 





v 

N 3 2 

6 A.s' + Bs! + ‘+ 
Sp orate Tae 


= ' ee t. ees: '. ' tevit. vi Vera ' _ We 
A Ni (m 4 eC Oe K ae cs Ks N vt ‘ Ys K *( 1. N 3 
r E 1g 


' 2 a ee ne = ' tevr.e _ ' ' _ wt. _ Tt - ie 
+ Ys N (IY K _) N. K oe 3 Ks (m x ant I. N <) 
rE 1g 1 
ee ' eo , ae = 3s ' t JS 7 1 vitae: ict Iw, 
2 N. (m a6 oe N. 1s K = + Ke N op Ks Nyy 5 
r 18 1g r 
' onl’ a ' '¢ ' ats = 1 feet 2a ' ite t. 
Ss Y¢ K Ne Y¢ Kas a N 2 Ys N oon Ys NY (Ty K ) 
r c 1g 1g 
Es ' tev _ ' wt, ' t _ owt, ' ' tf y Ve 
No 5 Ns Nt => a ‘y ee a5 7 a N 5) 
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- N, (m’ - YK, + Ng YK) + Ky ows + Ki, NOY’ + Y; BN 


- 6. VP v ¢ 6. vp > 
‘pnt oc vt weet ovit wert cnt plrevt owt pve t - vV're\N? 
3 KMD N Yi a oe Ys ae Ni iy 4 Ni Kp +K go Y Ny 
Key int. 
oe Vp 


' vty t wlve * pen? © vt wevt Cnt plvt — pt vine 
Ng ivNe + Rg ys ¥ Fs RNG Yves Ne vie Ky vg | 


CZ3436 
C 


Mw 


QAQAQ 
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APPENDIX C 


VERTLIN. FORT 


REAL A,B,C,D,E 
ENC eG eh vie Gay eee 


READCS, *)) WI XUD, ZWD, REY. Abd 
READC4Y,*) RNUD,Z2D,XWD,RMOD, RMWD 
READCY,*) X@D,ZUD,RMO,ZW,KkU,Z2 
READC4,*) RU, AW, XO, ZU, PME 
READCY,*) ZTHETA, ATHETA, RMW, ADE, ZDE, Rive 
REAL AiC5),A203),A304),A4 C4) 
COMPLEX = ZIC5)42262)5:23503) 224020 
UR Be 055-15) 
FORMATC IX, "IN WHAT FORMAT DO YOU WANT THE ANSWER 2 3 
DENS NSO AL WRITE 1',/7, 
: NONDIMENSIONAL URET ETO) 
READ C75 20) 209 
eG 0 aay Zl 
WRITECS,4) 
WRITE(6,4) 
FORMATC 1X, "OUTPUT VALUES ARE IN DIMENSIONAL Onn © 
GOr-TO..3 
WRT Ee S5o) 
URL TE Coes) 
FORMATC IX, ‘OUTPUT VALUES ARE IN NONDIMENSEONALTS OF 
CORTINUE 


POPMAT CII) 


XLU=XKL/U 
XLUZ=XLU*FXLU 
ALU3Z=XLU2*XLU 
XLUY4Y=XLU3S*XLU 
KLUS=XLUY*XLU 


A= (WI-XUD) *(WI-ZWD) *(RIY-RMOQD) -XWD* RMUD*ZEDSZUD*RiWD ae 
~(WI-ZWD)*RMUD*XOD-ZUD* KXWD* CRIY-RMNOD) SRNND* (WI SxU0? ae 


A ea ert Ree meme oer ea oe Loo oe eee Vo 
B= ~(WI-XUD)*(WI-ZWD )*RNC-Z2ZW* (NIT =kUD) ee 
—-XU*(WT-ZWD)* (CRIY-RMNOD) -XND*PMUD* (20). MUXXMD*ZOD 
—-RMUD*XW* ZOD-ZUD* RMND *XO-ZUD* RMW*ZOD-ZU~RMWD*XOD-( WT -Z2npee 
RHUD*XO-RMU* C(NT-ZWD) *XOD+RNUDAZW AK OR SAU ee 
~ZU*XWD* (RIY-RMGD)-ZUD*XW* (RIY-RMOD) -RMWD* (WI-AUD) *( 223 
=PMWE CNHIT-=XUD) “Z2ZeD + RUN DAU Aen 


C= -(WI-XUD) *(WI-ZWD ) *RMTHET+ZW* CWI-XUD) *RNO+AU* (WE Zee 
PMO+ZW*XU* CRIY-RMOD )-XND*RMUD*ZTHETA-RNU RD C2 
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VERTLIN.FORT (cont. ) 


RMUD*XW* (ZO+NT ) -ZUD*RMWD*XTHETA-ZUD* RMW*XO-ZUFRMWD*XQ-ZU*RMW* 
AeUS NT =ZUDJ*RNUD*XTHETA-RMU*(NI-ZND)*XStRMUD*ZWN*EXO+ 

2a eve UD eR che ecu AND Ret CUD XW* RNS 

ous ZU COREY “ROOD Ria (Wa xUD) * Ze t A-RMW* CWT S=XUD)* 

(ZO+WT )+RMWD*XU* (CZO+WNT) +XU* RNW*ZOD-RMU*KW*ZOD 


DO WNL Wh 


a 
D= Coen oe ON lawn) wu CWT GWD) RUTH EL=2W*xXU*RMS-RMNUX* 

MUD ZIRE TASRNUDAKH* ZIRETA-RMU PAN CZStWT J=ZUD*RMW*XTHETA 
acu ROD Tinta Kou Rh KOS RMU AC WI-ZWD) *XTHETA+RMNUD*ZN* 
AU wee RU et OUD rae PAs ZUD* XW RMIT ET*tXN*ZU*RMNG 
—RiW Ct =xUpD  ZIHETA+RNWDS XU Pt ZIRETAtRIWNtXU*t( ZO+WT ) 

623456789112 34567892 1234567803123456739412 345678951234 56789612345 

C 


Doe oon ou Ck hit aR AN ee a len et HET AtZW*RMUXXTHETA 
1 An Coen ti +XU* RAW ZIHETA 


C 
2345 


WD = 


ote co | 
Z| 


OQWitaHy 
Wot uo ey 
OQW0rmes 
% *%¥ * *¥G 
6 PS OG OMY 


QAQ 


Ze WRETE 


( 
¢ 
50 FORMAT OPFRPS.mOnr CHARACTEREISTIICGC ESUATION® -7) 
E 
paige AE Ah 2: 
51 FORMA pee CO le OD ale Kay ©?) 
pee 
se 
u 


100 PORT 


ro OQ tw 


Cc 


c 
CAneZPOELRCA],NDEG,Z1,IER) 
C23456789 | 


C 


WRIT 
200 FORM 


OO 
7~=OO 


TOOTS SOP PCHARACTERIS TIC ESUATION' 472 


—> ~~ ——b YO . 
>< BOBO MDD >< DY DO 


~ OO x OO-r 
eawwigwwy 


~ (hn 


Zr, 6 dF ORT REAL',11X%, ‘IMAGINARY’ ) 


1 
1 


Pm 

WW 

4 

ae | 
PM write ties 
ROW ROAWITOM 


E 

Zi 

= WRIT Z 

= 201 FORM 6. 
c 


ie) 
is 
2F i Snes) 


C 
See er CH RESPONSE TRANSFER FUNCTION NUMNERATOR POLYNOMIAL 


Te 


VERTLIN: FORT “(cone 


C 
C 
C 
ATHETA = RMDE*(WT-XUD)* (WT-ZWD)+ZDE*XWD* RNUD+XDE*ZUD*RNWD 
1 -RMNDE*XND*ZUD+KDE* (WI-ZND)*RMUD+ZDE iD Cul ee 
C 
C23 45697 
BTHETA = -RMDE*(WNT-XUD) *ZW-RNDE*XU* (CNT -ZU 00 7 eee 
1 ZDE*XW*RMNUD+XDE*ZUD*PMNN+XDE*ZU*RND-RHDE*XWD*ZU 
2 —-RMDE*XW*ZUD+KDE* (WI-ZWD) *RNU-X DE 2 Se 
3 =ZDE*RAUD*XUTZDE* ROW CSD 
C 
C 
CTHETA = RNDE*XU*ZW+ZDE* XW*RMNUtXDE*ZU*RUW RIDE 3 a 
1 Sk DE SZ PiU K=2 UE Seni 
C2345678 
C 
ERC S23 2435 1 
34 RIHETAS-ATHEDASX LUZ 
BTHETA =BTHETA*¥XLU 
C 
C 
32 WRITE CSS S52) 
WRG E 065.2) 
S2 FORMATC 1X; !’PLTICH POLYNONTAL COEFFICIENT Se 
Weir C5555) 
WEELEC 6539 
53 FORMATC7X>s “ATHETA' , OX, BIMETA COX. Clive Rae 
WRETECS,.100)5 AUWE DAs Bae eo beer 
WRETECG, 100) ATHETA -SIREIA, CINETA 
C 
NDEG=2 
AZC1)=ATHETA 
A2(C2)=BTHETA 
B2C3) =CTHETA 
C 
C 
CALL ZPOLRCAZ NDEG, 227 25R) 
C 
C 
WRETECS, 57) 
WRITEC6,54) 
54 PORMATUIA, ROOTS) 
WREEE CS 205) 
WELT EC6s 20a) 
WRITE CS 20s) 922 
WRETECG,.201) 822 
C 
~€ \VERTICAL VELOCITY TRANSFER FUNCTION NUMERATOR SPOu yhOur 
C 


AW = ZDE*(WT-XUD) *(RIY-RMOD)+XDE*RMUD*ZOD+RIDE*ZUD*XOD-=Z0E ee 


1 *XQOD+XDE*ZUD*( RIY-RMNOD) +NDE* (NT -XUD) Zep 
C234 5677 
C 
BW = -ZDE*(WI-XUD) *RMO@-ZDE*XU* (RLY-RMNOD) +FXDE*RMNU*2CD NDR ee 
= 1 *(ZO+WT)+MNDE*ZUD*XO+RNDE*ZU*XOD-ZDE* RUDE XC-ZDE oh 
= 2 —-XDE*ZUD*RMNG+XDE*ZU* (REIY-RMGD)+RADE “Chi Ue eee 
3 —~RINDE Au Ze) 


C2345678911234567892123456789312345678941234567895123456789612345 
CW = -ZDE*(CWIT-XUD) *RMTHET+ZDE*XU*RNO+XDE ©“RUUD*“ZIRE DA ee 


72 


VERE b bhi OR iam ccomt. ) 


] SCZ rien oe AUD nthe At hinDe ZU xZO-ZDE*RMUD*XTHETA 
oo on OK DES ZUDARPMIBET-“ADE*ZUt*RMO+RMNDE*CWT-XUD) 
Pee CT he PA SRM DE*XU*(ZO+tWT ) 
Cc 
C 
Dia ove su Ru lhno en DE <RMUSZTHETAPRIDE*ZUCXTHETA 
1 ~COon ene nth etna De ZUKRITRELTSORRDE“xXU*ZIRETA 
C 
Cc 
C 
IF(CJ) 42,42,41 
4 1 AW=AW*XLU3 
BW=BW*XLU2 
CW=CWFXXLU 
Cc 
Cc 
U2 WRETECS, 105) 
BrEeErEe C6, 1052 
iS Ment Cine VERTICAL VELOCITY COEFFICIENTS", 7) 
WRITEC5, 106) 
Viet e (6, 106) 
106 MORMATOC7/ A, AW',9X,"BW',9X,"*CW*,9X,  DW') 
WRITECS,100) AW,BW,CW,DW 
WRITEC6,100) AW, BW,CW,DW 
C 
C 
NDEG=3 
A3ZC1)=AW 
A302) =BW 
A303) =CW 
A3C4)=DW 
C 
Cc 
CiboechrOLR (AS, NDEG, 73; TER) 
Cc 
C 
Here to, 202) 
Me C6, 202) 
202 BORAT CI, “ROOTS* ) 
Beka EGon.2 05) 
WReet es (65205 ) 
Wiehe oa, 201) Z3 
Mee CO7 201) 23 
C 
C 
moo456789 
Dien KW > 2zWb) = C(CRIYSRNODI+RNDE*XWD*ZED+t+ZDE*RMND*XOD+RMNDE* 
—_ 1 tae oO D+ZDE*XWD*(REY-RHSDI-XDE*RMNWD*ZED 
Cc 
Cc 
Bg ae ADE CNT -ZWD)*RNS-ADE*( RIY-RMOD) *ZW+RMDE* XWD*( ZO+WT ) 
1 +RMDE*XW*XZOD+ZDE*RMWD*XO+ZDE*RMWXXOD+RMDE* CNI-ZWD) *XO 
Meee ce ZW *XOD-ZDE*XWD* RMNS+EZDE*KW* CREY-RMNOD)-XDE*RMWD 
: oe (2OFWT )—-XDE*RMIW* ZOD 
C 
C 
eee Oe CWT ZWD RMT ET ER DE*ZWNeRMOERMNDE*KXWND*ZTHETA+RMDE*XN* 
1 GZe Ln =ZDETRUWND* ATHETASZDE*RUN*AOTRMDE*CNT-ZWD)*XTHETA 
Z eran AO -CDE*KWD*RHITHET—ZDE*KW*RMNO-XDE*RMNWD*ZTHETA 


Q 


C2 


QQQ 


QQ 


QQ aa 


61 


62 
60 


70 


3 
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VERTLIN .PORT (eons) 


—XDEARHN* (Zo Sie 


DU = XDE*ZN*RMTHET+RMNDE*XW*ZTHETATZDE “Rin Aloe 
- RMDE*ZW*XTHETA-ZDE*XW*RNTHET-N DE Pa Zee 


PCI G2 2 Oi 
AU*XLU3 
BU*XLU2 
Cus 


oo 
aa 
Wout 


2tik>’ BU, IX s Cue sigi xe Oe 
U,BU,CU,DU 
U,BU,CU,DU 


WRIT 
WRIT 


we - 


Wo ul i ov 


PPPrsy 
FEE LYO 
oN ON ON ON DY 
FWNM 
Vw aviv ey |i 
OQ tr 

ccca 


CALL ZPOLRCA4,NDEG,Z4,IER) 


NN 
££ 


FORWARD SPEED POLYNOMEAL COERP ICE ae 


wo OH & F 


QaQAQ 


2, 


HORLIN. FORT 


Rea Ns D>» Ds E 

Pit eGoar NDEG, EER 

READC4, *IWT,RIX,RIY,RIZ,XL,U 

RPEADCis +a Von, LPHIL,YVD,YPDZYRD,YDR, YP 
READ(4,*)REV,RKP,RKR, RKPHL,RKVD, RKPD, RKRD,RKDR 
»RNR,RNVD,RNRD,RNDR,RNPD,RNP,RNPHI 
303)3,A204),A4C4) 

Mea 2)522030,24(3} 


"IN WHAT FORMAT DO YOU WANT THE ANSWER?"*,/, 
WRITE 1° 37; 
SIONAL WRITE 0%) 


GUZOr wre Z 
Pa pia 
he) 

t 


—/ 


7 OULEUL VARUES (ARE TN SPIMENSIONAL FORM’) 


— 


i] 


~ WOT BrP eV MHXs ew 


"OUTPUT VALUES ARE IN NONDIMENSIONAL FORM’ } 


XLU2 = XLU*XLU 
eGUS = ALUZ*XLU 
XLUGY = XLU3*XLU 
XLUS = XLUG*XLU 


A=(WT-YVD)*CRIZ-RNRD)*CRIX-RKPD) + RNVD*(-RKRD) *YPD 
Perr RD*RKVD*(-RNPDJ-RNVD*YRD*(RIX-RKPD)-RKVD* (RIZ-RNRDI*XYPD 
yee —Y VD) *(-PENPD) *(-RERD) 


Peo Ow tr VDD CREZ-RNRDDS RRP = RNR*(CWT-YVDI*CRIX=RKPD) 
fev Ghee RNR D) *(REX-RRP OD) tRNVD* C-RERD) F¥YP-RNVD*RER*YPD 
aoe OS RkeD)*YPD -~YRD*RRVD*“RNP +YRD*RKV*(—-PNPD) 
Bower | kk TD C-RNPD) + RNVD*YRD*REKP + RNVD*(WIT-YR)*(CRIX-RKPD) 
7 BNVELRD*CRIA-REPD} = RKVD<(RIZ-RNRDI*YP +RKVD*RNR*YPD-C RIZ-RNRD) 
Piov~ YP D 
SON - vv Dy C-RERDJ ERNPO S CWI -71 VD) *RRR*C-RNPD} 
oe ov C>RRPDI*(-RNPD) 


SHOU FWN— 


oe ee ~ CW VDF RKP =(WUT-yYVD)*+*(REZ-RNRD) *RKPHE+YV*C(RIZ-RNRD)*RKP 
RN ty VAC REX REPD)+RNVD* (-RERD)J*YPHE-RNVD*RKR* YP +RNV*(-RKRD) 
orev Pike Y PD -YRD “Rvp eRNPHESSyYRD*RKV*RNP +(WT-YR}*RKVD*RNP 
eee CHRP DD ST eRNVD*YRD*PRPHE = RNVD*(WT-YR)*RKP+ 
RNV*YRD*RKP + RNV*#¥CWT-YRI*CRIX-RKPD) - RKVD*CRIZ-RNRD)*YPHI+ 
RKVD*RNR*¥ YP - CRIZ-RNRDI*RKV*FYP +RKV*RNR<YPD +C(WT-YVD)*C-RKRD} 
Seteeeerdiam sy YD *RER*RNP = YV*¥C->RERED?*RNP - RKR*YV*C-RNPD) 


Dine wh — 


QQ 


OQ 


OAQ0Q 


QQ 


QQ 


96 


HORLIN.FORT (cont.) 


D = RNR* (WT-YVD)I*RKPHI + YV*(RIZ-RNRD) ©REPRDSR 
1 -~RNVD*RKER*YPHTI + RNV*(-RKRD)*YPHI - RNV*RER*YP — YRD*Re eee 
2 +(WT-YR)*RKVD*RNPHI + (CWTI-YR)J*RKV* RNP -RNVD* CWT =7 2) eee 
3 RNV*YRD*RKPHI -RNV*CWT-YR)*RKP +REVD*RNP eae 
G - (RIZ-RNRD)*RKEV*¥YPHI + REV*¥RNR*¥YP = (CUT=YVD) “Ree eee 
5 YV*(-RKRD) *RNERE +OReR yy Ree 
E = -RNR*YV*RKPHI ~- RNV*RER*YPHI + CWS Re eee 
1 -RNV*(CWT-YR)J*RKPHI + RKV*RNR*YPHD + RRR ere 
TRO) 22,227 211 
21 A = A¥*XLUS 
B = B¥XLUY 
C= CALS 
D = D*¥XLU2 
Bt= AL 
ZL WRITECS; 50) 
WREDE Coy 20) 
50 FORMATCIX, "COEFFS. OF CHARACTERISTIC EeuU- me ae 
WR 2 CS yo) 
WRITE(6,51) 
25 FORMATC7X;,°A*,» 14X,°B',14X,'C', 14X,'D" 7 16k 
WRITE CO; 100) (AS, C.D, E 
WRITES, 10:0) SA7 38> CC, De 
10 FORMAT Clk, SE 14.057) 
NDEG = 4 
AlC1) = A 
AAG 2) = 8 
AGS), = -C 
AlC4) = D 
APCS Jecs3. 5 
CALL ZPOLRCA IGHDEG,Z1;>IERD 
WRITE(5,200) 
WRITECO, 200) 
200 FORMATC1X;,*"ROOTS OF CHARACTERISTIC EOUlI 2G Gy) 
WRITECS 7205) 
WRITE GO, 2C 5) 
205 FORMAT CIOX, (REAL; 11%, "DMAGTINARY 
WRETECS) 201) 21 
WRIETEC6, 20 12) 
or FORMAT CIA AF lo. Gn 7) 
ABETA = YDR*(RIZ-RNRD)*(RIX-RKPD) - BRNDRMC> Ree ee 
1 = RKDR*YRD*(-RNPD) + RNDR*YRD*( RIX Seo 
Z + RKDR*(CRIZ-RNRDI*YPD -YDR*(—-RKRD) Clee 


BBETA = -YDR*(RIZ—-RNRD)*RKP —- YDR*RNR SCR pe 

1 -RNDR*(-RKRDJ*YP + RNDR*RER*YPD + REDE Dee 

2 + REDE*(WT-YR)*(=RNPD) SRNDR UR De Rie 

3 ~- (WT-YR)*RNDR*(RIXZ-RKPD?® + RKDR*( RES ee 

4 —- RKEDR*RNR*YPD + YDR*(-RKRD)J*RNP + DR eR Clee 


TC 


HORLIN.FORT (cont.) 


CEE PA = DR t CRIZ-RNEDIARKPHE + YDR*RNR*RKP 

1 -RNDR*(-RERD)J*YPHI +RNDR*RKR*YP + RKDR*¥YRD*RNPHI 

Zak R URS CW T=1R)*RNP See NDR*YYVRD*tRKPRo + CWI-YR) *RNDR*RKP 

See De eae) een ROR OR ee t COR C REED tRNPAL 


HD =1DR*RER*RNP 
623456789 
C 
C 
DBETA = YDR*RNR*RKPHI + RNDR*RKR*YPHI - RKEDR*¥*CWT-YR)*RNPHI 
ime Cee him nPHE SRNRAREDR*YPHIT = YDR*RKER*RNPHI 
ie CU, 32, 3 
al ABETA = ABETA*XLUG*U 
BBETA = BBETA¥FXLU3*U 
CBETAW— CBETA*XLUZ*U 
DBETAW= DEETA*KLU*U 
32 WRETECo7 a2) 
WRETECG, 52) 
aA FORMieMG@i Sie DESLEP “POLYNOMIAL COEFFS..(".7 3} 
URITE@s, 53) 
WRITE CO, 53) 
5 PORMapGrn, ABETA™, OX, BERETA’, 8X, "CBETA’, 8X, DBETA ) 
WRITECS, 10O0)ABETA,BBETA,CBETA,DBETA 
WRITEC6, 100)ABETA,BBETA,CBETA,DBETA 
NDEG = 3 
Ae C1) = RBEETA 
A2(2) = BBETA 
A2(C3) = CBETA 
A2C4U) = DBETA 
GCAnE ZPOLCR(CAZ, NDEGs ZZ, TER} 
WREEE CSO 4 } 
WRITEC6, 54) 
54 RORMAT CIA ROOTS.” 2) 
Wire 5.2 05) 
BW Reb i Gals 2 OS) 
Wee ieCs, 201022 
Mie Pee 2.0 122 
& 
C 
C 
APHI = RKDR*¥C(WT-YVDJ*CRIZ-RNRD) - YDR*RNVD*(-RKRD) 
ieteenw oe ReVD*YRD — RKDR*RNVD*YRD + YDR*RKVD*CRIZ—RNRD) 
Zee DR* (WI-yVDJ*C=RERD). © 
C 
C 
C 
Bemis = -RKDR*(WI-YVD)J*RNR = RKDR*YV*CRIZ—RNRD) 
ft ae NVD*RKR = YDR*RNV*(-RKERD) -— RBNDR*RKEVD*CWNT-YR) 
Zt DR Riv eR Det REDR*RNVD* (WE=YR) = RKDR*RNV*YRD 
Seek ny De RNR tt YDR*REV*CREZ“RNRDD + RNDR*CWT-YVDI*RKR 
P tote RNDIR~YV © C—-RKRD) 
C 
C 
Sence— PRR DR*YV*RNR + XYDR*RNV*RKR = RNDR*RKVECNT-YR) 
Pete Ra DR*RNVECWI-YR) = YDR*RKV*RNR =RNDR*YV*RKR 
IFOCJJ42,42,41 
4 APHI = APHI*XLU3 


Bene eB PHEtXLUZ 
CrPHr =-CPHI*xXLuU 


C 
42 


QAQ 


QQ 


RQAQQ QQQ 


QAA 


405 


106 


ZOZ2 


61 


bic 
60 


70 
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HORLIN.FORT (cont.) 


WRETE.C55. 105) 
WRITE Gon.103) 
FORMATC 1X, "ROLL POLYNOMDRL COEFFS = fe 
WRITE CS 100M 
WRITE CG, t0lc) 
FORMATC7X,*APHI’, 11X,' BPH’, )txXs) Ceri 
WRITE(CS, 100 )APHI,BPHI,CPHI 
WRITEC6, (O00 }APHT  SBPHIT, Crear 
NDEG = 2 
A3C1) = APHI 
A302)” =sS PRL 
ASQ 300 = (CruT 
CALL ZPOLR(A3,NDEG,Z3,IER) 
WRETE Co, 202) 
WRETE CG 5.202) 
FORTUALTGIX,. ROOTS») 
Weil oGs, 205) 
WET BCG. 203) 
WRETE CS, 2071923 
WRETE CG, 2013273 


APSI = RNDR*CWT-YVD)*CRIX-RKPD) + RKDR*¥RBNVD*YPD 
=——YOR*REVD*(=ENPD) + YOR*RNVD* CREA see 
= RNDR*=RAVD*YPD — RKEDRFCWTHYVD) + Ck ieee 


BPSI = ~-RNDR*(WT-YVD)*RRP - RNDR*YV*+ CRIX—Rir DD 
+ REDR*RNVD*YP +RRDR*ARNVAYPD +YDR*RE VO ie 
—~YDR*RKEV*C-RNPD) - YDR*RNVD*RERP +YDR*RNY Ce oe ee 
= RNDR*REVD*FYP = RNDR*RRVAYPD + REDR* CUT 
+ REDR*YV *C—“RNPD) 


SI = -RNDR*(WI-YVD)*RKPHI + RNDR*YV*RKP +REDR*RNV DY Poe 
RKDR*¥RNV*YP +¥YDR*¥RKVD*RNPHI +YDR*RKAV“RNP = YDR*RNVD Ree 
YDR*RNVEREP = RNDR*¥REVD*YPHE = RNUR 

REDR* (CWI-Y VD) *RNPHI = RRDRALV ARNE 


DR*YV*RKPHI + RKRDR*¥RNV*YPHE + YOR REV RNeoe 
*RKPHI - RNDR*RKV*YPHT “SRE “2 See pee 


"APSI' »11X,°BPSI',11X,'CPSi’ pal a ee 
JAPSLT, BPSE,Crsh, Vest 


Je) 


HORLIN.FORT (cont.) 


PRG COs OO JAP Se, Bhool, Chou, DPS © 
NDEGe=— 3 

AGC1) = APSI 

AY(C2)} = BPSI 

Ao Cs — CPS. 

AGYC4) = DPSI 

CALL ZPOLRCAY,NDEG,Z4,IER) 
WIRE GE 5, 202) 

Whe 6720.2) 

WRG 55200 

Were 6, 205) 
Were CS 2012 
WRITEC6,201)24 

SOP 
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APPENDED 


AXIAL FORCE 


: 2 2 : : = 
m| i - vr t+ wq - x, (q +r ) + yo (pq r) + 20 (pr +a) | 


| p Af 2 t 2 7 ! | 
fe ee > aa +X r“ + X r 


ro 3 . ¥ t t ] 
tS ry [x.' a + a wg 


ae 4?| X ree Xt ves 'w? | 
2 uu vv ww 


D 2 2 r t 2 t 2 - 1 z 
ioe see LX su gr’ 887 +X oq gn" 88" +X op op 6b? | 
ten Sale 2 
+3 pt ja, +b, uu. tc, a. ] 
- (W - B) sin @ 


gee ee § 2 4? 


2 
“ K oy brn z 


Pp. | if t 
_ 2 ss aT ww) m 


2 2 
PT oCen Ca 1c | (m1) 


a9) 


LATERAL FORCE 


m| v - wp tur - yo (r? + p?) + 26 (ar - p) +x, (ap +#) | = 


tho fyo rey pty = ge re r | 
3 r DF ES otgen) 2 PY Soy EN repo 


t t 
ars 


eh Ply vey ay Cae 
ee Mie v vq we 


eee jot 
o> ee Cet, UP tt 


| . ot 
ice ujr[é6r + aac [(v* + w)2 lirt | 


- i 
+= L° pes u” + Dean UV Sal Vv I(v? +w?)? | | 


28 2 { 1 1 4 ] 
+> 4 ee. vw t ¥, tue br 


+(W - B) cos 6 sin 6 


p 


wary ee oe Pur (9-1) 


n 


Bere | 2)2 Tne 
ts x L vn rane te iG ee oe a | (1-1) 
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NORMAL FORCE 


m| w - uq + vp “te (p? +97) +X, (rp - V+YVe (rq +9) | 


Pp 4 | : 2 1 pe ' | 
se i Comet) aa p a r + 2,5 se 


Feed fra (pags as 'vr+Z ' vp | 

2. Ww Vl vp 

ea eZ, 'ulq|és +°Z LS Nee 2y2 

2 q 2 “1a1ss wip 
; Iw l 








: 2 
L° Es u? + Zi uw Daren, w|(v" + w?)? |] 


1 
‘ £° L210 ulw| + are Iw (v? + wf | 


pos [ 2 2. 2 | 
5 L a v~ +Z u~ 6s + Zen u*~ -6b 


és 


+ (W - B) cos Ecos 


P_ 43 : 
+= ¢ Zon ug (7)-1) 


ae 
wlo 


: b§ 
2 1 1 2 22 a ah = 
L | an uw t Z win wi(v? + w*)? | 4 Zeon O<% |n-a) 


10) 


ROLLING MOMENT 


pe (I - 4) qar+ (r + pq) LS + ? aia) le + (pr - q) Nose 
tml yo (w -uqtvp) - Zc (v - wp tur) |= 
+t sl K. p+K,'i4 ee Cee Pat Kp Piel 
+> “lk Pp +K_'ur + Key | 


'p af j ' 4K t r | 
es L ev vg US wp one 


ome 
ag 7 Ky! a? +K fav + Ky! v |(v? +w?)?| | 


. 3 [ i 4.2 ] 
ao L jK VM + K u% 6r 
Bs (Yo Wee Yp B) cos @cos $ - (2. W - z,8) cos § sin > | 


P ,2 
+2 + Ky, u® (N-1) 
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PITCHING MOMENT 


I ¢ I -] olen oo ae nese 
yo + G1.) tp - (pt ar) 1p 


tm [2g (a- vr + wa) - x0 (Ww - ug + vp) | = 
p il Pa fae 1 2 1 2 1 l 
ty ELM) a+, Pet My!) + MP 4M ala | 
fae iM. w4M lovee t NO Ny ] 
2 Ww vr vp P 
p Mi 2 2\5 ] 
+— M 1 2 
5 £ ; Sa eG lies ulalés + Mii oily + w*) lq 


} 
8) 3 bse 4 \ 2 2,2 
to 2 |My! u? + ML ow +My [(v? + w?) | | 


w | 


i 
hPL Mitte] eM lw tv? +0 |] 


1 | i aye mie 
i £ MY +M,ou 6s + M,. 


ee OW - xp B) cos 9 cos } - (z. W - z, B) sin 8 


G 


0D 4 { 
+= 447M u -1) 
5 an q (n 


Se 2, = ; :) 
a : + § -! 
3 LM Oe allen wli(v we ) l+ Moen au (7-1) 


O5 


YAWING MOMENT 


; : _¢ ye ee 
\ 7 oe 1) Pq - (q + rp) Sy Dr eee? Aiea Pp) i 
A - ° =) 
+m | X>, (Vv - wp + ur) - y¢ (u - vr + wq) | = 
eae [Not +Ni'p +N 'pq +N__'ar ] 
ee tae p P+Nog' PatNa,' a + Nyy, | ‘Flt 
2 1 ays \ i ' ] 
eS £ N° fs) ae OE Gas whe vq 
soft s*[N | ap tyNe are WN 'ulr|ér +N ' |(v? + w2)? | ] 
Z Pp es lr | 6r [vir : 
1 
p ,3 2 2 2,2 
ta- £ [Nt u +N Ca ee l(v? + w | | 
+= “21 N 'vw+N.! u br | 
vw ér 
+ (X, W - x, B) cos 8 sin ¢ +(y. W -y, B) sin 6 
o t4N_,! ur (0-1) 


p ,3 2\2 2 
+ r : n - 
5 | N pe sl v|(v* + w*) | +N, bu? {(n 1) 


<. 
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KINEMATIC RELATIONS 


U? =u? + vv +w* 


-u sin @+vcos @sing¢+wcos @cos ¢ 


Ne 
it 


¢ =p+y% sin 6 


9 _q- cos @ sin ¢ 
cos ¢ 


r+ 6 sin ¢$ 
cos @cos¢ 
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(edasgn SOL a7 
LSYT Ss! LIA‘ Awl L /WALSAS/ 
CEYNOON /SLNANT/Z 


-(C9°O)ANTIWS 
“(AL VOTSOxd cto C21 O1SX be Cer Sock (2) xi a0 n 37 


MMNH’QMNH’MANH‘SAANH‘ AANH 
“YANHSOANH‘AANH‘S LOGANH ‘ANH’ SNH ‘4 MANH S OMNH SYOYNH 
“LOGUNH‘YNH‘ MYONH ‘DdNH‘! LOGGNH‘ ANH‘ AYGNH’YGNH /AABIOON/ 


AMMWHSMMWH’OMWH ‘ AMWH ‘ LOOMWHY 
“MYMWHSMWH‘S AAWH 4 YMAWH ‘AAWH‘SWH ‘ MUWH ‘dd uWH ‘OUOWH 
“FOWH’ SGOWH 4 LOQOWH ‘OWH 4‘ ddWH‘ ASGWH‘SUGWH’MYWH /ASSNOW/ 


UMYH’ UMHS MAMH ‘AAYHS OAH‘ LOGAWH’ ANH‘ SSH’ SHH 
“LOGUYH ‘YNH ‘MOH’ OdWH’ ddWH’ LOGAWH‘ dWH‘HQYH /AATOINS 


IMMZ4MMZSOMZ‘SSAMZ‘S LOOMZMYMZSMZ‘SAAZ‘SUNZS AAZ 
soe ouzZ diz 307° SU0Z" LOGO7 Uz ddzZ ASdzZ SdzZ av Zz 7 ISA00Z7 


UMKS AMR 
“MARS AAAKSAAKSUARSOARSSARS LOGARSAXSSHA ‘SAMAK 
“HYCUASTLOGHAS HAS HORS OdK  ddk’ LOGdK‘ dk‘ SHAK‘ MAK /AATHOORK/ 


AMMX MAX OMYSAAAXSAAXSHAXSANX’! LOGNYX 
“Hux’dadux’ vox’ asasax’sdasax’ ayuduax’/ywaudx‘ Id‘ lv /aagQOX/ 


XZIA‘ZAKIG’XXIA 
‘ZIG’ RIG’ XIG‘ AZ‘ AR‘ AX’ 9OZ‘ OR OX AG‘ TA‘ NA /SSYW/ 


NOWWOO 


NOWWOD 
ey 
NOWWOQD) 
6BB8LISPECZO 
Lf 
NOWWOO 


NOWWOD 


NOWWOO 


NOWWOO 


NOWWOO 


NOWWOD 


NOWWOOD 


NOWWAO 
68L9OSPECO 


(Z-O’H-V)B8a1V9N LIOITdWI 
qOWwndL ANTLAOXENS 


mM o+ CNX) a aero y 


“SOa INIMOTION SHEL YOd CTSIOW YVANTTINON 300-9) 
HLAYL SHL SO SAATLIVATHAG ALVIS S3LVINOTVO ANILNOYHNS SIHL 
ROO Gok kkk ok oktokgkok dk ktokokdok kkk ktok ok & 


- 


d XIQGNHddV 





TAgonw 


mr OU O 


nd 
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LOGUYH*ASTZH = XZId= = (9 FIWL 


AXI1S=- = (S‘b)WL 
LOGdYH*SG1ZH |= XI = Cb‘ H)WI 
QkeWA = (E‘°P)WL 
LOGAWHRKEIZY = OZeWHe = (27° DI) WL 
LOGOZ*B1ZN = OX*WHe = (S’°EIWL 
Qk*eWG = (bE) WT 
LOGMZ*+E€1ZH =- WH = CFELE)WL 
LOGNA*E TZN = OX¥WA = (947) WI 
LOGdA*¥DIZH = OZeWH- = (b6‘'Z)WOL 
LOGAK*¥€7ZN |= WA = (2‘Z)WI 
DheWwh- = (9%) WL 
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